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INITIALIZATION, CONCEPTUALIZATION, AND APPLICATION 

IN THE 

GENERALIZED FRACTIONAL CALCULUS 


Carl F. Lorenzo 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

Tom T. Hartley 
University of Akron 
Department of Electrical Engineering 
Akron, Ohio 


Abstract 

This paper provides a formalized basis for initialization in the fractional calculus. The intent 
is to make the fractional calculus readily accessible to engineering and the sciences. A modified 
set of definitions for the fractional calculus is provided which formally include the effects of 
initialization. Conceptualizations of fractional derivatives and integrals are shown. Physical 
examples of the basic elements from electronics are presented along with examples from 
dynamics, material science, viscoelasticity, filtering, instrumentation, and electrochemistry to 
indicate the broad application of the theory and to demonstrate the use of the mathematics. The 
fundamental criteria for a generalized calculus established by Ross (1974) are shown to hold for 
the generalized fractional calculus under appropriate conditions. A new generalized form for the 
Laplace transform of the generalized differintegral is derived. The concept of a variable structure 
(order) differintegral is presented along with initial efforts toward meaningful definitions. 


1. Introduction 

The development of fractional calculus is nearly as old as the development of the actual 
calculus. Discussion as early as 1695 by Leibnitz, and later by others, indicates an interest in 
integrals and derivatives of fractional (non-integer) order. Excellent historical summaries of the 
development of the fractional calculus can be found in Oldham and Spanier (1974) and Miller 
and Ross (1993). Although much attention has been directed toward the development of the 
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fractional calculus, this has largely been in the domain of the mathematician and the applied 
mathematician with relatively little work being done in the engineering and applied sciences. 

Strong motivation exists for the study, development, and generalization of the fractional 
calculus. This may be readily validated by looking ahead to the Applications section of this paper 
where a broad range of problems are described and solved using a generalization of the fractional 
calculus. The unifying concepts and notation of the fractional calculus provide a significant 
benefit that greatly simplifies the solution of certain partial differential equations (distributed 
systems). Perhaps the strongest motivation to develop the fractional calculus is the belief that a 
wide variety of physical problems and engineering disciplines that have resisted compact (and 
first principles) description and solution when using the integer order calculus will yield to the 
methods of the fractional calculus. These include such areas as heat transfer, boundary layer 
behavior, and many problems where the major recourse has been to probabilistic methods. 

This paper looks at some of the primary problems with the fractional calculus, as it is now 
embodied, and attempts to rectify some of these with a modified (alternative) embodiment. The 
basic approach that has been taken in this work has been to make the defined mathematics as 
maximally applicable to the problems of engineering and science as possible. To this end, basic 
distributed dynamic systems have been developed for use as reference fractional systems. The 
authors have tried to maintain maximal generality in these reference systems (Hartley and 
Lorenzo (1998)) as required for the new definitions. 

The problems that the authors perceive to bar widespread application of the fractional calculus 
in the engineering sciences will now be discussed. The operative basis definition for this section 
will be the contemporary, Riemann-Liouville definition, which is that integration of arbitrary order 
is given by 

where -°°< c <t and v > 0 and T(v) is the gamma function. The notation here is as follows: the 
left term is that of Ross (1974) and, originally, Davis (1936), and the middle term is that of Oldham 
and Spanier (1974). For simplicity in this paper v is constrained to be a real number (or real 
variable). The reader is cautioned that the notation used later in the paper is the same as Ross but 
has a different (generalized) meaning, which is defined there. 

Contemporary differentiation is defined as 

iff) 

,o;/(0=-^,o, ■’/('). (i.2) 

where v = m-q and v > 0 and m is an integer > v . For the reader new to the fractional calculus 
it should be noted that several other basis definitions have been forwarded by various authors. 
Reference to these may be found in the following excellent publications: Miller and Ross (1993), 
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Oldham and Spanier (1974), Samko et al. (1993), and elsewhere. It is further noted that equation 
(1.1) is a positive running ( t increasing) integration since it is required that t- r > 0. 

Alternatives are discussed in the literature that may be viewed as negative running integrations 
such as the Weyl definition for fractional integration (see for example. Miller and Ross (1993)). 

One of the fundamental problems of contemporary fractional calculus, is the requirement that 
the function f(t ) and its derivatives be identically equal to zero for t = c (i.e., the time of 
initialization) (Ross (1974)) or lacking this, to limit the functions handled to special classes (see 
for example. Miller and Ross, (1993)). This is needed to assure that composition (the index law) 
holds, that is, to assure that 

lD " , £>,"/(')= cDrfit )■ d-3) 


It is difficult, in the engineering sciences, to always require that the functions and derivatives 
be zero at initialization. This fundamentally says that there can either be no initialization or 
composition is lost. Thus, it is not in general true that 

d Q d Q f 


f- 


df Q dt 


Q 


= 0 


(1.4) 


(see for example, Oldham and Spanier ( 1974), pp. 155 and 82-87). 

Thus, when solving a fractional ordinary differential equation of the form 


( 1 . 5 ) 


additional terms must be addended to equation ( 1 .4) 

/ '~~o - C,t Q *‘+Cit Q ' 2 +. . . (1.6) 

dt v d v 

to achieve the most general solution 

/ = +C,,c- I +C ! t <! ! + . . . c., 1 ( 1 . 7 ) 

dt v 

The reader is referred to Oldham and Spanier (1974) for a detailed exposition in this area. The 
added terms are analogous to the effects of the constants of integration in the integer order 
calculus. 

Both of the above issues are related to the inadequacy of the composition law. This 
inadequacy of the composition law, in fact, points to a problem with the fundamental definitions 
used as the basis of the fractional calculus. 
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Flowing from the above inadequacies, the Laplace transform of the fractional differintegral, 
currently (Oldham and Spanier (1974)) given as 





s k f (0), 


df 


( 1 . 8 ) 


is found to lack sufficient generality in the initialization term(s) for many applications. 

The objective of this paper is to provide a formalized basis for the initialization of the 
fractional (generalized) calculus so that it may be readily accessible to the sciences and 
engineering. In particular, a set of definitions will be provided which formally (explicitly) include 
the initialization functions. Conceptualizations of fractional derivatives and integrals will be 
shown. Applications and physical examples of the basic elements (many drawn from electronics) 
will be presented along with examples from dynamics, material science, and more to indicate the 
broad application of the theory. Also, the fundamental criteria for a fractional or generalized 
calculus established by Ross (1974) will be investigated. Those criteria which apply to the 
proposed approach will be proven and rationale for compromising/generalizing those remaining 
will be given. The Laplace transform of the generalized differintegral will be derived. Finally, 
some physical conceptualizations will be given for fractional differintegration and some 
applications will be used to demonstrate the mathematics. 

The authors have attempted to write to both the mathematics and the engineering/scientific 
fields. Thus, the engineering reader will find familiar block diagrams associated with the 
mathematics and, hopefully, the mathematician will find sufficient rigor in the proofs. 


2. Motivation for Initialization 

In real applications, it is usually the case that the problem to be solved is in some way 
isolated from the past. That is, it should not be necessary to retreat to — °° in time to start the 
analysis. Usually, the analyst desires to start the analysis (for instance, in the case of a system of 
ordinary differential equations (ODE's)) at some time, say t 0 , with knowledge (or assumption) of 
all values of the function and its derivatives. Specifically, /(r 0 ),/ , (r 0 ),/"(r 0 ),.../ ( ” ) (r 0 ),in the 
case of ODE’s. In modem parlance, this collection of constants is called the system state and 
contains the effect of all the past history. 

The desire to initialize fractional ordinary differential equations (FODE’s) or a system of 
FODE’s continues to exist in the fractional or generalized calculus domain. To understand the 
application of fractional calculus in terms of fractional or extraordinary differential equations, the 


NASA/TP— 1 998-2084 1 5 


4 



authors considered a system (Hartley and Lorenzo (1998)) which can be described as a semi- 
differential equation or semi-integral equation. Physically this system is the semi-infinite lossy 


line, which is described by a diffusion equation of the form 

= a 3 = v (0,f) = V/(f). v(~,0 = 0, V(x,0) given. (2.1) 

dt dx r dx 


Here, v is the voltage, i is the current, V,{t) is a time-dependent input variable, and a is a 
constant, which is one over the product of resistance per unit length, r , and capacitance per unit 
length, c. 

One way the behavior of the semi-infinite line can be described is in terms of its input 


behavior (impedance) at the open end of the line, that is, as a semi-differential equation. 

However, to practically use such an FODE to describe the behavior at the input, requires the 
addition of a function of time. In terms of the physics, this time function relates back to the initial 
voltage distribution (distributed initialization) on the semi-infinite lossy line. The details of this 
derivation are given in Hartley and Lorenzo (1998). The primary result for the voltage at the 
input of the lossy line is given by 


V(0.s) 


r/(0,s) 


s 

a, — 


\e\4 V(A,0) dX , 

o 


( 2 . 2 ) 


where A is a dummy variable of integration and s is the Laplace transform parameter. 

Equation (2.2) was determined using a conventional iterated Laplace transform approach 
applied to the diffusion equation. However, an attraction of the fractional calculus is the ability to 
express the behavior of the line (a distributed system or mathematically partial differential 
equation) as part of the system of distributed equations using fractional ordinary differential 
equations (FODE’s). Such a fractional (extra)ordinary differential equation for this application 


would be of the form 


d m v(t) 


dt 


1/2 


'4a i(t). 


(2.3) 


assuming v(*,0) = 0. To initialize this distributed system, a function of time, l f/(t), must be 
added to account for the integral term in equation (2.2). Thus, 

^-+y f(t) = r4a i(t). 


d l/2 v t 


dt 


1/2 


(2.4) 


The focus of interest is the integral term or yr(t) of equation (2.4). Clearly, one can addend 
such terms in an ad hoc way to the FODE’s which are being solved. However, it is the authors’ 
opinion that a better approach is to formally include an initialization function, y/(t ) , in the 
definition of the fractional differintegrals. This is the approach that will be developed later in the 


N AS AjTP— 1998-208415 


5 



paper. Clearly, if the analyst is constrained that the initial function value and all of its derivatives 
are zero, the range of applicability for this entire class of problems (which will probably include, 
eventually, all distributed systems) will be greatly limited. It is the authors’ contention that all 
fractional ordinary differintegral equations require initialization terms to be associated with each 
fractional differintegral term, in order to complete the description. Indeed it will be shown that 
some ordinary differential equations (i.e., integer order differintegrals) may require a 
y/(t) initialization. This requirement is a generalization to the requirement of a set of 
initialization constants (the state) in ordinary differential equations. Fundamentally it is the 
information required to start the integration process of the differential equations while properly 
accounting for the effects of the past. Before presenting the formal definition of this generalized 
calculus, it is useful to examine some interpretations of fractional differintegrals. 


3. Conceptualization 

The interpretations of fractional differintegrals presented in this section are based on the 
contemporary definitions see, Ross (1974) and Oldham and Spanier (1974). The extensions of 
these interpretations to the definitions that follow in this paper will be obvious. 

3.1 Riemann-Liouville Differintegral 


3.1.1 Interpretation 1 


The Riemann-Liouville definition for the fractional integral is 

r,/\i -, /(0-=nJ(»-Tr7(T)rfr <?>0. 


This definition is extended to fractional derivatives by simply differentiating; thus 


(T2_ f 

[d(i-o)r 



q>0, m>q. 


(3.1.1.1) 

(3.1. 1.2) 


and m is an integer. 

A graphical display of the semi-derivative and semi-integral of f(t) = cos(0.4;rr)is presented 
in figure 3-1-1. For this case a - 0, and equation (3. 1.1.1) becomes 


d U2 f{t) 1 r cos(0.4ttt ) 

[*wr = rfi) dt 


(3.1.13) 
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Plate (a) of figure 3-1-1 shows the fixed part of the integrand; this is unchanged for all semi- 
integrals regardless of f{t), t is a parameter on this plot. 

Plate (b) is the time function to be semi-integrated, plotted as a function of T . The complete 
integrands, which are the product of the functions in plates (a) and (b) are as shown in plate (c); 
again various values of t are shown as parameters. The semi-integral, plate (d), is the locus of the 
values of the integrals of all parametric values of t as shown in plate (c), (and more). The semi- 
derivative is obtained simply by differentiating the half integral, and is shown in plate (e). 



tau 



Figure 3-1-1. — Interpretation of fractional differintegrals based on Reimann-Liouville definition. 
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3.1.2 Interpretation 2 


The Riemann-Liouville integral can also be viewed as a convolution integral as 


d-fjt) 

K-or 





[ f f(*)dr 


q > 0. 


( 11 . 2 . 1 ) 


In this instance both f{t) and h{t) = 1 / r(q)t~‘ l+ ' are considered to be causal (i.e., no convolution 
response can be obtained before f(t) is applied), and * represents the convolution operation (see 
for example, Churchill (1958)). An excellent discussion of this interpretation of the convolution 
integral is given in Gabel and Roberts (1973). This viewpoint is reflected graphically in figure 3- 
1-2 where again /(/) = cos(0.47rt) and q = 0.5. Here, plate (a) shows h(z) versus T , the negative 
function h(- r) is shown in plate (c), and h(t- r) for t = 5 is plotted on plate (d). The function 
/(r), to be fractionally integrated, is presented in plate (b). Plate (e) shows the full integrand, 
h{t- t)/(t) for t = 5.0 . Finally, the integral of this product becomes the t = 5.0 value of the 
fractional integral, plate (f). 



Figure 3-1-2. — Riemann-Liouville integral interpreted as a convolution integral. 


NASA/TP— 1998-208415 


8 





3.2 Griinwald Differintegral 


3.2.1 Time Delay /Conveyor Interpretation 


Another, more physically based, interpretation of a fractional differintegral can be made 
based on the Griinwald definition. This definition, considered as the primary definition by 
Oldham and Spanier (1974), is given as 


t — a 


d*f(l ) 
[d(l-a)]' 


GRUN 


= lim 

N— >«x> 


N ) 

r(-«) 


AM 

I 


j=0 


r U-<t) 

r(y + l) 




(3.2.1. l) 


where q<0 indicates fractional integration and q > 0 indicates fractional differentiation. 
Understanding of this interpretation will be increased by studying a particular case. Consider the 
much simplified approximation to equation (3.2.1 .1) for q = 1/2, N = 4. Then from equation 
(3.2.1. 1) 


[d{t-a)T 



r(- i / 2) 

r(i) 


f(t)+ 


r(i / 2) 
r(2) 

' r(3/2) 
r(3) 


f\t- 


4 )' 


f 



r(5/2) 

r(4) 




(3.2.1.2) 


Each of the function evaluations is a time delayed evaluation of f(t) . Further, each of these time 
delayed functions is multiplied by a constant, which is r(j-q)/ r(-q)r(j + l)and also multiplied 



In diagram form, it is seen that the base function f(t) is evaluated unshifted and 
progressively time-shifted in increments of [(t — a) / N) for ( N — l) many times. In this case the 
diagram is as shown in figure 3-2-1. Thus, the semi-derivative of /(/) is seen to be a summation 
of progressively delayed evaluations of f(t) multiplied by progressively decreasing constants 
and finally multiplied by cr' l/2 . Now N = 4 is a very crude approximation to the semi-derivative 
and equation (3.2. 1 . 1 ) calls for the limit as N approaches infinity. It is instructive to allow N to 
approach infinity. 
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Figure 3-2-1. — Time delay approximation to semi-derivative. 


The following case is considered, <7 = 1 / 2 , N = 10 , 000 . Then equation (3.2. 1. 1 ) is much 
more closely approximated by 

r(-i/2) ] 

r(i) + 
r(i/ 2 ) ( ( t-a X\ 

r(2) '(* uooooJJ + 

r( 3 / 2 ) ( J t-a ^ 

r( 3 ) J [ t uooooJJ + 


J t-a 
f t — ~h 

•'l UooooJJ 


J t-a 

f t-2 + 

uooooJJ 


d m f{t) 1 10000 J r ( j 1 / 2) f ( t-g V) 

[d(t-a)] U2 ~ r M] r(; + l) f V y ll0000jj + 


T(9997.5) ( ( t-a X) 

T(9999) ^l ?- 9998 llOOOo]J + 


T(9998.5) ( (t-a 

— f\ f_90QQ 

r(ioooo) [ l 10000 


It can be seen by comparison with the previous approximation, that the gamma function 
based coefficients are the same for the first N terms. The time shift factor (incremental delay) 
({t-a)/ N) is of course very much smaller and indeed approaches zero as N goes to infinity. 
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A diagram for this approximation is given in figure 3-2-2. The diagrams of figures 3-2- 1 and 
3-2-2 are easily associated to a physical interpretation by considering the behavior of a conveyor 
figure 3-2-3(a). The analogy considers a conveyor with some material (e.g., sand) dropped from 
a hopper to the conveyor surface. The height of the sand on the conveyor is related to /(/) . 
Height sensors are evenly spaced such that the time for the conveyor to move the distance 
between two sensors is a ={t- a) I N when the conveyor moves at speed V(/). The figure shows 

a time function f{t ) and its shifted counterparts f(t-o),f(t-2o) f(t -(N - l)cr). The 

outputs of the height sensors are multiplied by appropriate gamma function coefficients as in 
figures 3-2-1 and 3-2-2 and summed, etc. The result of the summation then is the fractional 
differintegral. This analogy evenly distributes N sensors between x = 0 and x = t , therefore the 
spatial spacing changes as t increases. 

An alternate conveyor analogy allows a fixed spatial distribution if in the Griinwald equation 
(3.2. 1.1) the substitution N =(/- a) /AT is made. In the limit as AT ->0, AT is always chosen 
so that N is an integer, then this gives 


DH{t) = lim 

« * J ' ' A T* v/1 


(AT) 


-H N - 1 


AT— >0 




r(-«) ji rO+i) 


j(AT)). 


(3.2.1.4) 


This situation is indicated in figure 3-2-3(b). Assuming no initialization, the sensors at 
fU - mo) at time t = 0 will show no output. It is clear that preloading the conveyor to the right 
of the zero point is the same as initializing the differintegral and is the same as extending the 
function backward in time. 



Figure 3-2-2. — High order approximation of the semi-derivative. 
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Figure 3-2-3(a). — Conveyer analogy: Fixed N, here n = 1000, sensors uniformly spread 
between 0 and t, spacing varies with t. 
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1. J /(/) — |o !': 2 ' S2 Uniformly Spaced Sensors 

y i*^*s — 

I I I I I I I I I I I I I I I I I I I 1 

V „ = 1 ft / second 








Figure 3-2-3(b). — Conveyer analogy: N variable in t, fixed in space, 
sensors uniformly spread in space. 
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A third conveyor analogy is possible. In this case the weighting sensors are evenly distributed 
and fixed in space. For this case the velocity V(t) of the conveyor is not constant. It varies from 
infinite speed at the start (/ = 0), to zero speed at t = °° , in such a manner as to satisfy the 
required time delays. (No figure shown.) 


3.2.2 Geometric Interpretation of Fractional Differintegral 


A geometric conceptualization may be obtained based on the Griinwald definition. Namely, 


„£>,»/(/)= lim 


teTV r 0-«) 


r(-«)£S r 0 + i) 


/ 


t-j 


t-a 

N 


(3.2.2. 1) 


For simplicity we will write {t-a)/ N = AT , then since N — > °°, as AT — > 0, 




(3.22.2) 


The nature of the definition may be explored by considering the adjacent j th and 

(j + l)th terms. In a general sense if the terms are additive and q < 0, then an integration is being 

effected. If the terms are differenced and q > 0 then differentiation may be suggested. Then 


U 


£>,"/(') = lim 

' J ' ' Ar->o 


h 


r(j-q) f{t-jAT) 

T{-q)r{j + 1) AT q 


+ 


r(y+i-g) /Mt+Qat ) 

r(-q)r(j + 2) AT q 


+••• 


\. ( 3.22.3 ) 


Dividing through by the coefficient of the j th term, and combining the yth and (y + l)th terms 
gives 


•ww- 


. _ 1 /V 

f f(*-j&r)+P /(^-0+ 1 ) A7 ’) 


rCc 

AT‘< 

H r 

/ 


(3.22.4) 


where 


Tjj+l-q) Tjj + l) Tjj-q) 

r (j-q) r(y+ 2 ) an a r(- 9 )r0+i)’ 


j = 1,2,». , 


and a is the coefficient of the j th term. Plots of a and (3 are presented in figure 3-2-4. Now 
using the gamma function relationships this may be simplified to /3 = [j -q)/{j + l) when 
j — q > 0 and j + 1 > 0. 

We now restrict our attention to 1 > q > -l.Then the simplification for fi substituted into 
equation (3.2. 2.4) holds for all terms with j > 1 . Further, under the q restriction, (3 is constrained 
to the range of l>(3>0 and (3 = 1 when q = -1 and approaches 1 as j — > oo for the remaining 
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values of q (see fig. 3-2-4). Thus, it is seen that so long as the terms above j = 1 exist, (3 will 
always be positive; hence the signs of the adjacent terms will be the same. Therefore, all of the 
considered cases will involve integration processes! It is noted that the sign of a only influences 
the sign of the resulting value of the integration. 



order, q 



Figure 3-2-4. — Plots of a = T (j - q)l T (- q) T (j + 1) 
for various j and /} = (/- 9 VO + *) versus j. 

Now equations (3. 2.2.2) and (3.2.2.4) are applied to consider four special cases. 

Case 1: q = -1 

Now q - - 1=>J0 = 1, and a = 1 ; further, equation (3. 2.2.4) holds for all j. Equation 
(3.2.2.4) then simplifies to 

a D;' fit) = fimj- • •+A7’[/(r - jAT) + f(t - (j + l)Aj)]+- • • } (3.2.25) 

which is readily seen to be two terms of a conventional integration. 
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Case 2 : q = 1 

When q = 1 => a = 0 for all j except = 0 and 1 . From equation (3.2. 2. 2), the j = 0 and 
j = 1 coefficients are determined to be 1 and -1 respectively, giving 


' ,J AT^o\ A T J 


(3.2.2.6) 


which is the definition of the conventional derivative. 


Case 3 : 0 > <7 > — 1 

For simplicity we will start with q = -/ 2 => fi_y 2 = (y >2) / (7 + 1) ; thus, for all 
j > 0 => 1 > P _y > 0 . Here (3 has been subscripted to indicate the value of q . It is also true that 
1 > a > 0 for all q in the range (fig. 3-2-4). Thus, 

u D~ y ' f(t) = Jim {• • -+a A T*[f(t-jAT) + p_ yi /(r-(y + l)A7)]+ •••}, (32.2.1) 

and for a general q in the range 

„ 0,V(<) = i™,[ "+« AT— '[/(/ - jAT)+P t 1 )AT)]+ -•} (3.2.2.8) 

Thus, since both a and (3 are always positive, this summation will be seen to be an integration 
process (fractional). A geometric approximation to this integration is shown in figure 3-2-5(a). 
Consider q = - 1 , then ^ a f{t - AT ) A T is an area represented by that area beneath the 
a f(t - jAT ) curve. If q is taken as - 2 (outside the domain of consideration here), then 
^ a f(t - AT)AT 2 would be a volume. Then the series /(/ - AT)AT q for 0 > q > -1 may be 
thought of as a “fractional area” (or “generalized dimension”). 



Figure 3-2-5(a). — Geometric interpretation of fractional integration; integral is 
area under 0 (fi t - jAT) curve multiplied by AT 11 ' 1 . 
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Case 4: 1 > q > 0 

Again, for simplicity, we will start with q = X and generalize. For q = }< => ^ = ( j - Yi)l (j + 1) , 
here for all j > 1 => 1 > ^ > 0 . It can be seen that will be positive for all when j > 1 . It is also 
true that 0 > a > -1 for all q in the range (see fig. 3-2-4). 

Thus, 

„ D'{- f ( t ) = Jim {• • -+a A T ~' /: [/(r - jAT ) + Py 2 f(t -(j + 0 A T )]+' ‘ * }♦ (3.2.2.9) 

or in general for q in the range 

u D//(/)= lim {•••+« AT- <l [f(t-jAT)+^ /(t-(;' + l)Ar)]+ - }. (3.2.2.10) 


So, after the first j = 0 term, it is seen that all terms again are a direct sum of negatively 


weighted functions, again an integration process. However, the effectiveness of the weighting 
A7” v is changed since now q > 0. The first j = 0 term for all q is AT~‘ l f{t)\ thus, considering 
the first two terms (eq. (3. 2. 2. 2)) yields 


„D, */(>)= Hni 




Ar->o 


AT q 


(3.2.2.11) 


This brings in an effective differentiation (for q > 0 ) though scaled by AT q instead of AT as in 
the case of order 1 differentiation. Therefore, fractional differentiation may be viewed as a 
combination of integration and differentiation. 

If <7 is taken as 1 in equation (3.2. 2. 2), then, of course, the equation returns the rate of change 
of the function. Taking a value of q = 2 for the exponent of AT (outside the range of 
consideration), then equation (3. 2.2. 2) yields precisely d~ f / dr , an acceleration. Thus, for 
1 > q > 0 the terms of equation ( 3 . 2 . 2. 1 1 ) can be considered as a “fractional rate” of change of the 
function. Figure 3-2-5(b) shows the j = 0 and j = 1 points of a geometric approximation to the 
qth fractional derivative. The slope between the curves multiplied by AT~ q+] is, loosely, a 
geometric interpretation for this part of the fractional derivative or “fractional rate , the 
remaining terms are interpreted as in figure 3-2-5(a). Of course, this interpretation becomes 
meaningless in the limit as AT — > 0. 



Figure 3-2-5(b). — Geometric interpretation of fractional differentiation,) = 0 => j = 1 . 
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3.2.3 Geometric Interpretation 


Considerable insight may be obtained from the examination of adjacent terms of the 
Griinwald definition of the fractional differintegral. The geometry of both fractional integration 
and differentiation at the infinitesimal level are contained in the form 


a 


f(l-jAT)+Pf(t-(j+\)AT) 


(3.2.3.1) 


The graphical interpretations are as shown in figure 3-2-5. 


4. Definitions for the Generalized Calculus 
With Initialization Functions 

In view of the motivations presented in the earlier sections of this report and the desire to 
apply the fractional calculus to many practical applications, the following definitions are offered. 
Discussion will begin with a consideration of the Riemann-Liouville integral. In the development 
that follows, attention is restricted to real values of the order, q, of the various differintegrals. 

4.1 Riemann-Liouville Basis 

Two types of initialization are considered. In the first type it is assumed that the differintegral 
operator can only be initialized ("charged") by effectively differintegrating prior to the "start" 
time, t = c . This will be called terminal charging (terminal initialization) by analogy to the 
driving of a distributed semi-infinite electrical line, at its open terminal, by a driving function 
fit), for ( a<t <c) . This is contrasted to side charging (side initialization), the second type of 
initialization, where by similar analogy a fully arbitrary initialization may be applied to the 
differintegral operator at time t = c . Side charging is analogous to setting an arbitrary (voltage or 
current) distribution on a distributed semi-infinite electrical line. 

Terminal charging is considered first. It is assumed that the fractional integration of interest 
"starts" at t = c (i.e., point of initialization). Further, /(/) = 0 for all / < a , and the fractional 
integration takes place for t>c>a. 

The standard (contemporary) definition of a fractional integral will be accepted only when the 
differintegrand f(t) = 0 for all t < a . 

Then, 

„ D: fit) = /(*) dr, V > 0 , t>a, (4.1.1) 
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subject to /(f) = 0 for all t < a. The following definition of fractional integration will apply 
generally (i.e., at any t > c ): 

e D?f{t) = 77 -T }(f-r ) v ' f(z)dz+ y/(f-v, a, c,t), v > 0, t>c, (4.1.2) 

r(v ) \ 

c> a and f(t) =0 V t< a. 

The function y/(f -v,a,c,t) will be called the initialization function and will be chosen such 
that 

u Dl q f(t)= for t > c. (4.1.3) 

Mathematically, then, for t > c and v > 0, 

a } (f - t)' - ' f(z)dz = c D7 f(t ) = } (t - r)’"' f(r)dt + y{f -v,a,c,t). (4.1.4) 

t c t 

Clearly, since J f{z)dx = J f(z)dz + J f{z)dz , 

a a c 

\l/(f-v,a,c,t) = a D'f(t) = p^-yj(r-T)‘ - 7(TVT t>c, v>0. (4.1.5) 

This expression for 1 j/{t) applies only for the terminal charging condition. Clearly, y/ brings to 
the definition of the fractional integral the effect of the past, namely, the effect of fractionally 
integrating f{t) from a to c. This effect will, of course, influence behavior after the time, t = c . 
The y/ function has the effect of allowing the function f{t) and its derivatives to start at a value 
other than zero, namely, the value 

u D: q f{t) L, 

and continues to contribute to the differintegral response after t = c. That is, a function of time is 
added to the uninitialized integral, not just a constant at t = c. 

The integer order integrals under terminal charging are of special interest. Evaluating 
equation 4.1.5, for example, for v = 1 , indicates that y/(f ,-l,a,c,f)=constant. The general case 
is readily shown to be 

n 

y/(f,-n,a,c,t) = ^ l c i t' , n = l,2,3--- . 

/=0 

This, of course, is the same effect as seen in the integer order calculus using the “constant of 
integration.” 

When side charging is in effect as opposed to terminal charging, then, 

y/ = y/[t) (i.e., is arbitrary). 
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It is important to note that the initialization of the qth fractional integral of f{t) is not unique in 
the following sense. That is, /(/)can be considered as a composite function, for example, 
f(t) = g(t)U(t-c) + h(t)(U(t-a)-U{t-c)), where U(t) is the unit step function 

JO r<0 

[1 t> 0. 

Then for this composite function f(t) , it is the function g(t)U{t - c ) that is being 
differintegrated and h(t)(U (t - a) - U {t - c)) is the function on which the initialization is based. 
This is analogous to choosing an arbitrary constant value to initialize (the integration of) dy / dt 
in the solution of an ordinary differential equation. The point is that in many (perhaps most) 
applications it will be the analyst's choice as to how h(t)(u(t- a) -U{t-c)) will be chosen in the 
initialization of t . Df f(t) . In many cases, for example for the integer order (differintegrals) 
calculus as considered in the contemporary context, h(t)= constant will be the appropriate 
choice. Issues of discontinuities at t = c of course must be addressed. 

To extend the definition to the fractional derivatives, some issues must be addressed. The 
definition of the fractional derivative raises the following important questions in the context of 
initialization. Do fractional derivatives require an initializing function in general ? Further, do 
integer order derivatives in this context require initialization functions? Clearly, as we 
commonly think of derivatives, in the integer order calculus, the derivative is a local property and 
is represented geometrically as the slope of the function being considered and as such it requires 
no initialization. In the solution of differential equations the initialization constants which set the 
initial values of the derivatives really have the effect of accounting for the integration of the 
derivative from minus infinity to the starting time of the integration (of the differential equation). 

The authors’ study of the semi-infinite bar dynamics (Hartley and Lorenzo (1998)) together 
with implication from the conveyor analogy based on the Griinwald definition of the fractional 
derivative, indicate that the fractional derivative is not a local property as appears to be the case 
for integer order derivatives (in the integer order calculus). Further, if it is desired to use the 
fractional calculus (and hence the fractional derivative) to write FODE’s such as that representing 
the input behavior of the semi-infinite bar, then an initialization function is required to handle the 
effect of the distributed initialization (initial strain distribution in the bar). The input 
characteristic (for displacement to force) in this case is df(t)/ dt (that is, an integer derivative 
but in a fractional context). Therefore, in the opinion of the authors, the answer to both questions 
is clearly, YES. Physical examples will be shown later in this paper to further justify this 
position. The impact of this is to require an initialization function for the fractional (and integer 
order) derivatives. 
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Thus, a generalized integer order differentiation is defined as 

» m 

cDTf(t) = —f(t)+y/(f,m,a,c,t) t>c, (4.1.6) 

at 

where nt is a positive integer and where m,a, c,r)is an initialization function. This is, of 
course, a generalization of the definition of the derivative, and for many cases, for example, 
usually in the integer order calculus, i/r will be taken to be zero. It will be shown in a later section 
that, for m= 1 with the condition of terminal charging, y/(f ,l,a,c,t) = 0 . 

Now the uninitialized generalized (fractional) derivative is defined as 

a D*f(t)= y D? a D; p f(t) q>0, t>a, and f(t)=0 \ft<a (4.1.7) 

and (for convenience) m is the least integer greater than < 7 , and q = m — p. Now as in the 
fractional integral case y/(f ,-p,a,a,t) = 0. Further, since y/ (h,m, a, a, t)= 0, where 
h(t)= i D; p f(t) , this definition specializes to the contemporary definition of the fractional 
derivative. 

Now the initialized generalized (fractional) derivative is defined as 

,Dlf{t)= <D'," C DT f{t) q > 0, t>c>a , (4.1.8) 


where (for convenience) m is the least positive integer greater than q and q = m - p . 

The initialization of the fractional derivative under conditions of terminal charging will 
proceed in a manner similar to the fractional integral, that is, it shall be required that 

£>/ f(t)= U D? f(t) \/t>c with< 7 > 0 . (4.1.9) 


Specifically, this requires compatibility of the derivatives starting at t = a and t = c, for t > c. 
Then it follows that 

, D "‘ ( D; p f(t)= (1 D; p f(t) , q> o, t>c>a. (4.1.10) 

Expanding the generalized integral terms. 




• 1 

77 -T f ( t - t) /,_i f{x)dx + \i/(f,-p,a,c,t) = 

r lP) J 

\ ' "1 

— f (t - )dx + v (/ -p,a,a,t) 

r {p)i J 


t > c . (4.1.11) 


Since y/(f-p,a,a,t) = 0 and using the definition (eq. (4.1.6)) for the integer order derivative, 
equation (4. 1 . 1 1 ) may be written as 

r r [ (/ - t )'’" 1 f{x)dx+ y/(f-p,a,c,t)\ + ,m,a,c,t)= 

ApK J 

I t 

—7^\(t-xy- X f{x)dx +y/(h-,,m,a,a,t ) , 

r {P)t 


dt' 1 


d m 

dt" 


t >c. 


(4.1.12) 
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where /i, = a D t ''/(/) and h 2 = u D, p f(t). The integer derivative is uninitialized at t = a , therefore 
\ff(h 2 ,m,a,a,t ) = 0. Then combining integrals gives 

d m \ 1 c ^ 

i ]/(h ] ,m,a,cj) = — — r-T\{t-T) p ~' f{T)dT-y/(f-p,a,c,t) t>c. (4.1.13) 

V M^/ a y 

Under the condition of terminal charging of the fractional integral (eq. (4.1.5)), the argument of 
the derivative above is zero thus y/{h ] ,m,a,c,t ) = Oand it is seen that integer differentiation 
cannot be initialized through the terminal (terminal charged). This feature will also be 
demonstrated in the applications section to follow. For the case of side charging y/(f p,a,c,t) 
is arbitrary. Thus it can be seen from the above equation that either 
y/(f-p,a,c,t) or y{h\,m,a,c,t) can be arbitrary but not both while still satisfying the 
requirements of the initialization (eq. (4. 1.1. 3)). 

Thus, the generalized (fractional) derivative, side charging case, can now be stated as 

«-AV(0= , ' f(^)^ + ¥(f -p,fl,c,/)|, q>0, t > c, (4.1.14) 

and m is the least positive integer >q with q = m— p, or equivalently as 

(jin If d™ 

. D 1 /(0= ~J t „, J(f - t) p ' f(r)dr +-j^y/(f ,-p,a,c,t)+y/(h,m,a,c,t) q> 0 (4.1.14a) 
where m is as above, t > c,and h(t) = <( D“ /J /(r) . 

In the case of terminal charging of the fractional integral part of equation (4.1.14) 
y/(f,-p,a,c,t ) will be as defined by equation (4.1.5). 


4.1.1 Simple Example 


A simple example will be helpful. Consider the semi-integral of f(t) = {t - a)U(t - a ) , then 
for the uninitialized semi-integral (Oldham and Spanier (1974) pp. 63-64) 

u D; l,2 (t-a)U(t-a) = “ 0.7522 5(/-a) 3/2 , t > a. (4.1.15) 


Now, initializing from the point t = c 

I i 

t D, U2 (t-a)U(t-a)= Y^j^j(t-T) ]/2 (r-a)dr+li/(f-l/2,a,c > t) , t > c, (4.1.16) 


2 

t D) n {t-a)U{t-a) = ^ ((t-c) U2 (2t + c-3a))+v(f-\/2,a,c,t) , t > c. (4.1.17) 
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Consider now the terminally charged initialization, 

if/(f r W2,a,c,t) = ^^yJ(t-T)" 1/2 (r-a)r/T, t > c, (4.1.18) 

\lf(fr\/2MX\t)=^^—j[(t-c)' l2 (-2t + 3a-c)-(t-a) m (-2t + 2a)\, t > c. (4.1.19) 

For specific numerical values, for example, a = -1, c = 1 , the above results give 

A D~ u \t-a)U{t-a)= 0.75225(r + lf\ t>-l, 

? / x x (4.1.20) 

, D r{t-a)U{t-a) = ^=((/ - 1)' /2 (2f + 4)) + yf(f -Ml,- 1,1, t ), t > 1, 

with 

=^=(-{t-l)' n (2t + 4)+{t + \) U2 (2t + 2)), t > 1. (4.1.21) 

The numerical evaluations of these quantities are shown in the graphs of figure 4-1-1. It is 

observed, at least for this case, that 

limy/ (f,-v,a,c,t)=0. (4.1.22) 

f — 

Further, for t>c= 1 , y f{t) appears to decrease monotonically. 



Figure 4-1-1. — Sample problem: semi-integral of (t-a) versus Time, with a = -l, c= 1. 
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Summary of Definitions 


The generalized (fractional) integral, for v arbitrary, real, and non-negative: 
cD,'f(t) = -pjr-j J'(t-f)' ' f{T)dT+\j/{f-v,a,c,t) , v > 0 

t >c > a and /(/) = 0 t< a . 


(4.1.23) 


The generalized integer derivative: 

ini 

,DT f(t) =—f(t)+\ff{fm,a,c,t), t>c, 

at (4.1.24) 

where m is a positive integer and where y/(fm,a,c,t) is arbitrary. 

The generalized (fractional ) derivative, q and p real: 

i m i t j in 

< D? = £>;", D; p f(t)= — r(p) J(/-ry’-7(rVT +—\ff(f,p,a,c,t)+\f/(h,m,a,c,t), (4.1.25) 
where t >c, h- a D~ p f(t ) , q = m - p>0, p> 0 and m is integer. 


General case side charging : 


Wit) is arbitrary. 


Terminal charging case for generalized (fractional) integral: 

y/{f-v,a,c,t) = u D]'f(t)= -J— J(r-r) vH /(T)</r t>c,v> 0 


(4.1.26) 


for integer derivative 


V(h, m,a,c, 0=o. 


(4.1.27) 


4.2 Grtinwald Basis 


This discussion parallels that for the Riemann-Liouville basis. Again the "starting" time for 
the differintegration is t = c . An initialization is introduced to account for previous history and 
goes back to t = a , with f(t) = 0 for all t < a . Then for arbitrary q , 


,D;'f(t) = 


<i q f(t) 

[d{t-a)] q 


_ v - — - V f { t { t ~ a 

~ r(-<?) JrO + i) 7 ! J [ n 


t > a, (4.2.1) 


subject to f(t) = 0 for all t < a . 
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The following definition of Griinwald differintegration will apply generally (i.e., at any t > c > a): 

d ^ 1 ‘ ¥a{. fq>a,c,t) (4.2.2) 


cDlf(t) = 


[d{t-c)\ 


q I GRVN 


and f{t) = 0, V t < a , c > a. 

Again i f/ G (f ,q,a,c,t) is chosen such that e D l 9 f(t) will produce the same result as a D? f{t) for t>c. 
Mathematically, 


,D'/f(t) = c d?M 


+ W if q,a>c,t)= u D?f{t) , 


(4.2.3) 


GRUN 


for all t > c,and f{t) = 0 for all t < o,andc > a. 

Therefore, 

VW(/ ^a,c,t) = a D? /(/)- c d 1f{t\ GRUN , 

subject to f(t ) = 0 V/<fl and a < c <t 
or 

Vgrun (/ ,<? ,a,c,t)= a d‘f f(t)\ GRUN ~ c d, f(ty GRUN * 
subject to/(t)=0 V t < a a<c<t. 

Now substituting the Griinwald series definition from equation (4.2.1) gives 


(4.2.4) 


(4.2.5) 


y/ G ( f,q,a,c,t) = lim 


( \~ q 

t-a 


v N, j 




r{-q) fi r(j+i) 


^ t — a^ 


lim-! 

Af,-i 


(t-cV 


kNi) 2 ^ 4 ^- j\ 


f t-c ' 


r (-<?) m r(y+i) 


l (v 2 


(4.2.6) 


subject to t > c, and f(t) = 0 for all t < a . 

The above expression can be simplified greatly by making the delay increments in the two 
summations equal, that is, AG = ((t - c) / {t - a)) AT,. Then after considerable algebra and 
manipulation of indices 

<4 ' 2j) 

where AT =(t-a)/N ,, t>c>aar\dN\ and Ni are integers such that Ni={(c-a)/(t-a))Nu 
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4.3 Notation 


It is important to clarify notation before proceeding further. The following notation will be 
used in the remainder of this paper: 

c D?f(t)= c d?fit)+y/(f,q,a,c,t), t>c, (4.3.1) 


where c D)‘f(t) represents the initialized q th order differintegration of f(t). 
Thus, 


d‘/f(t) = 


d q fit) 


[d{t-c)]‘ 


= | 


ROSS 


OSlS 


(4.3.2) 


is an uninitialized generalized (or fractional) q th order differintegral and is the same as defined 
in Oldham and Spanier (1974) and Ross (1974). Its basis is either the Riemann-Liouville integral 
(and successive differentiations) or the Griinwald series definition. The notational form 

c D? f(t) , (4.3.3) 

is now used to represent the initialized differintegral, that is, it includes the y/(t) initialization 
function. It is unfortunate that this form is widely used in the literature at this time to represent 
the uninitialized differintegral; but it is the authors’ judgment that any confusion arising from 
this change will be more than offset by the convenience of the compact notation. As described 
above ys(f ,q,a,c,t) is an initialization for the q th order differintegral of /(f) where it is 
required that a < c and that /(/)= 0 for t <a. For convenience y/(f ,q,a,c , /) will occasionally 
be replaced by t f/(t) or y/(f,q,t ) where the meaning is clear from the context. In similar 
manner, d) 1 f(t) will occasionally be replaced by d q f(t). 


5. Criteria for a Generalized (Fractional) Calculus 

5.1 Ross’ Criteria 

This section examines selected properties, which can reasonably be expected to apply to 
fractional differintegral operators. A minimal set of criteria for a viable generalized calculus has 
been established by Ross (1974b). This set is accepted in principle by the authors and is quoted 
as follows: 

The Ross Criteria: 

“1. If f{z ) is an analytic function of the complex variable z, the 
derivative L D' z f(z) is an analytic function of v and z . 
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2. The operation c D' x f{ x) must produce the same result as ordinary differentiation 
when v is a positive integer. If v is a negative integer, say v = -n , then 

cDTf(x) must produce the same result as ordinary n-fold integration 
and c Dx f(x) must vanish along with its n - 1 derivatives at x = c . 

3. The operation of order zero leaves the function unchanged: 

///(A I - HU (5.1.1) 

4. The fractional operators must be linear: 

, D?[af(x)+bg{x)] = a < Df f(x)+ b < Df g(x). (5.1.2) 

5. The law of exponents for integration of arbitrary order holds 

cDx ,D7f(x)= cD;rf(x). ” (5.1.3) 

It is observed that the notation, u D‘‘ x fix), is exactly that of Ross and refers to the uninitialized 
differintegral. 

The requirement here is to examine the proposed generalized calculus against the spirit of 
these criteria. It should be noted that there is a minor conflict contained in these criteria. Namely, 
criterion 2 basically calls for backward compatibility with the integer order calculus while 
criterion 5 calls for the index law to hold. It should be noted that, for the normal Riemann- 

Stieltjes integration theory in the integer order calculus, 

d- m d m f{x)*d° f(x)=f{x) (5.1.4) 

for all f{x) and all m . This is observed directly by application of the “Fundamental theorem" of 
calculus (Apostol (1957) p. 215). This states that under appropriate conditions 

j/'(0*=/M-/M. (5J-5) 

C 

and it can be seen that the reversal of differentiation by integration really differs from fit) by 
/(c) , that is, by the initialization. Thus, the fact is that the integer order calculus itself does not 
satisfy all the criteria. The failure of the fundamental theorem to provide true reversibility of 
differentiation by integration is handled in the integer calculus in an ad hoc manner, specifically, 
by the use of the constant of integration and by the use of the complementary function in the 
solution of differential equations. 

The issue to be addressed is: Do the initialized versions of the Riemann-Liouville and Griinwald 
differintegral definitions satisfy the above criteria? Since Oldham and Spanier (1974) show the 
equivalence of the Riemann-Liouville and Griinwald uninitialized definitions, and since the 
Riemann-Liouville basis relates most closely to the classical calculus it will be the main 
consideration of this effort. Further, since the focus of this paper is limited to q in the real domain, 
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only the last four of the criteria will be considered here. It is noted that, at least for continuous 
functions, the condition f(t) = 0 for V t < a , at t = a , is equivalent to the condition that is used by 
Ross (1974a and 1974b), namely, f(t),f'(t),f"(t), ... f {n] (t) = 0 for Vn and Vt<a . 

A general consideration of the Ross criteria relative to the generalized calculus proposed here 
allows the following to be established. 

For ... f la) (t) = 0 forVn and Vf < a and under the condition of terminal 

charging, criteria 2,3, 4, and 5 hold for all D? f(t), with t > c. 

Proof: Since each criterion has been proven (Ross (1974a and 1974b)) for the uninitialized 
differintegral under the condition f{t),f'(t),f"{t),...f {a) (t) = 0 forVn and Vt<a.And 
since for terminal charging, the generalized integer derivative becomes 

■“/('). t>c, (5.1.6) 

at 

and further, because y/(f,q,a,a,t) = 0 Vgthen, 

=„£>;/(' toB (5.1.7) 

(Note „ D," /(/)| R0SS here represents the uninitialized differintegral). Hence, since all criteria hold 
f or „D?f(t)\ R0SS when (t) = 0 V n and V t < a , and since 

c Dff(t)= a D?f(t), for all t>c and all q, (5.1.8) 

it follows that the criteria also apply to the initialized differintegral, 

c Df f(t) , for all q and t > c. 


thus completing the proof. 

The criteria are now examined separately, to detail the proof for each criterion, to determine 
the implications relative to the initialization functions (y/ ’s), and to determine the limitations for 
the generalized case, that is, for the case of side charging and/or the case of an initialized integer 
derivative. 


5.2 Criterion 2: Backward Compatibility 


The purpose of Criterion 2 is to require backward compatibility with the integer order 
calculus. Namely, it is required that 

d n fit) 


Jf, fit) = 


df 


n > 0, 


or 


n = integer. 


(5.2.1) 


i: K-.~a »<<>. 
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Ross (1974b) uses induction to prove compatibility for the uninitialized fractional integral, 
thus when v becomes a positive integer n , 

r («)o 

produces the same result as n -iterated ordinary integration. Now by definition, 

cD\ n f(t) = :p^y l f(T)dT + y/{f-n,a,c,t) , n > 0, (4123) 

t >c>a, and f(t) = 0 V t< a, 

and since for side charging y/(f ,-n,a,c,t) is arbitrary, for backward compatibility (i.e., to satisfy 
eq. (5.2.1)) it may be taken as zero. For terminal charging, f(t ) can be defined to be identically 
zero for a < r < c . Thus, backward compatibility is obtainable for the generalized integral when 
v is a positive integer. 

The backward compatibility of the generalized derivative c D?f(t), q > 0 to the integer 

order derivative — — — , when q is replaced by n, requires that the zero property (criterion 3) 
dt" 

hold for the fractional integral part of equation (4.1.8). This property is independently established 
in the next subsection and will therefore be used freely here. Then as p — > 0, equation (4.1.8) 
becomes 

lim t D^ f{t)= lim ( D, m C D~ P f(t) 0, t>c>a 

q— 1 f>—> 0 

i m 

= c D?f{t) = —f(t) + il/{f,m,a,c,t) (5.2.2) 

and for terminal charging it has been shown that y/(f ,m,a,c,t) = 0 , and under side-charging 
conditions y/ is arbitrary and may be set to zero. Thus, backward compatibility for the fractional 
derivative is obtainable. 


5.3 Criterion 3: Zero Property 


For criterion 3, it is necessary to prove that 

cD? fit) = fit). 

Now, by definition, and using the Riemann-Liouville form for the fractional integral 


I 1 ' 

c D°f{t) = lim -pj-y \(t -z) y ~' f{r)dT + y/(f ,v,a,c,t) 


t > c. 


(5.3.1) 


(5.3.2) 
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and for terminal charging 


= lim 

v ->0 


-!-}(/■ r r 1 f(r)dz 9'"' Mdx 


r(v)J' - ' r(v) 

, A°/(')= lim^y j(r - zT' f{r)dT=\\m( a D]f(t)), 


and 


t > c. 


(5.3.3) 

(5.3.4) 


Now in this form the differintegral has the initialization, /(f) = 0,forV/ < a , that is, 
y/(f-v,a,a,t) = 0. Now Ross shows by several methods that 

j > 

Km rTtl ^ ~ ^ ' f^ dT =/W • 

v - >0 r(v ) j 0 


(5.3.5) 


Thus, it is seen that the initializing function does not affect order zero behavior of the fractional 
integral with terminal charging conditions. 

Approaching zero order from the derivative (positive) side, taking q = m- p=l- p , since 
no generality is lost by taking m= 1 , 

, . 0 , 7 ( 0 = lim , D, ’/(<)= lim, D\,D; " f(t) (5.3.6) 




.. d 

= lim — 

p-* i dt 


■p^J ( r_T ) p ' f( T )d?+v{f -P,a,c,t) +y/(h,m,a,c,t), t>c. (5.3.7) 
Then for terminal charging of the integer derivative y/(h, m) = 0 and for the fractional integral, 

j c 

= t>c, (5.3.8) 


thus. 


• D ° /(,) - ~ r) "" 

^~]f(r)dr. 


and 


t > a . 


(5.3.9) 

(5.3.10) 


Under conditions of continuity of f{t) and the order of operations can be interchanged 

and by the fundamental theorem of calculus 

t 

,.o,7(0= f/'Mdr =/(0-/(«)=/(0 . (5.3.i i) 

a 

since by definition f(a ) = 0, completing the proof. For side charging the zero order operation 
will return 

’ (5.3.12) 


, D?f(t) = f(t) + ¥ (t). 


That this is an acceptable and indeed important result will be demonstrated in a later section. 
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5.4 Criterion 4: Linearity 


Linearity becomes somewhat more involved with the presence of the initialization function in 
the definition. The requirements for the fractional integral will be determined first. 

The following two fractional integrals are defined: 

c D;' f(t) = f(r)dT+y/{f-v,a,c,t), t>c, and (5.4.1) 

r(v)J 

C D< g{t) = g(r)dr+if/(g,-v,a,c,t), t>c, (5.4.2) 

uvj, 

where/(f) = 0 Vf < a and g(t ) = 0 Vt < a. 

Now, for t > c consider, 

, D~' ( bf{t ) + A:g(/))=-p^y j (t - t)“' ( bf (r) + k g(r))dz + y/(bf + kg,-v,a,c, t ) 

=-^ \(t-rr' f(r)dr +-^-z \(i g(z)dr +y/{bf + kg-v,a,c,t ). (5.4.3) 

r(v) J T(v ) ' 

Now using the definitions from equations (5.4.1) and (5.4.2) to replace the integral terms in this 
equation gives: 

, D;'(bf(t)+kg{tj) - b{_ 

k( r D;’g(t)-ll>(g-v,a,c,t))+y/(bf + kg,-v,a,c,t), l>c. (5.4.4) 

Clearly, 

c D;'{bf(r) + kg(t)) = b [ D;'f(t) + k t D;'g(t), t>c, (5.4.5) 

and linearity is obtained if an only if 

i j/(bf + kg,-v,a,c,t) = b\l/(f,-v,a,c,t)+ky(g-v,a,c,t), t >c. (5.4.6) 

Fundamentally, equation (5.4.6) presents the requirement on the initialization of the individual 
functions to produce the same result as that given by the initialization of the combined functions. 

Extension of the linearity requirements to fractional derivatives, while not complex, is 
simplified by use of the following shorthand for the uninitialized fractional integral, 

f(t)= ( d; q f(t)=^\{t-Ty-' f(r)dr, q>0 (5.4.7) 
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and 


m> 0. (5.4.8) 

at 

The problem here is to show 

. D;(bf(t) + kg(t)) = b c D;f(t) + k c D; g (t) , r> 0, t>c. (5.4.9) 

Using the definition for the generalized derivative (eq. (4.1.8)), 

, D ; f(t)= ,D, m , D;' fit ) , t>c, (5.4.10) 

where m is an integer and m > v > 0 and r = m-v. Now, r D, m is an initialized integer order 

d"‘ 

mth derivative operator composed of two parts, the normal derivative operator, and 

y/(h,,m) the initializing function. Also, for equation (5.4.10), h l = a d~ r f{t). 

For t > c, linearity holds when for equation (5.4.9), LHS (left-hand side) = RHS (right-hand side), 
where 

LHS= D;{bf(t)+kg(t)). (5.4.11) 

By definition of the generalized derivative, with m and v as defined above, 

LHS = D m D - v (bf{t) + k g(/)) . (5.4.12) 

By definitions of D m and D ~' , 

LHS =d m (d-'(bf(t)+kg(tj)+\i/{bf + kg,-v))+y/(h,m) 

=d'"d-' {bf(t)+kg{t))+d m li/(bf + kg-v) + \i/{h,m ) , (5.4.13) 

where h(t)= a d;' (bf(t) + k g{t)). 

Using the linearity of the fractional integral and the ordinary derivative 

LHS =bd m d- v f(t)+kd m d- r g(t)+d m \if(bf{t)+kg(t)-v)+y/{h,m ) . (5.4.14) 

For the RHS of equation (5.4.9), 

RHS = b D r f(t) + k D r g{t) (5.4.15) 

and by definition of the generalized derivative 

RHS = bD m D~'fit) + kD m D~'git). (5.4.16) 

Now, by definitions of D m and D v and rearranging terms, 

RHS=bd m d-' fit)+kd m d-'git)+bd"'\ffif-v)+kd m \i/ig-v)+by/ip,m)+kii/il,m), (5.4.17) 
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where 


(5.4.1 la) 


p{t)= a d;'f(t), and l(t)= a d; v g(t). 

Therefore, linearity of the derivative requires LHS = RHS , then from equations (5.4. 1.4) and 
(5.4. 1.7), 

d m y/(bf(t)+ kg(t),-v)+\i/(h,m)= bd m y{f -v)+ k d"‘ii/(g-v)+ b\f/(p,m)+ k\f/(l,m) , (5.4.18) 

and requires satisfaction of equation (5.4.6), linearity of the fractional integral. Equation (5.4.6) 
can be differentiated to yield 

d'"y/(b f(t) + k g(t)-v) = bd"'ii/(f-v) + kd"'\f/(g-v), t >c. (5.4.19) 

This result substituted into equation (5.4. 18) gives the simpler requirement 

y/(h,m) = by/(p,m) + ky/(l,m) t>c , (5.4.20) 

for the fractional derivative. 

The results (eq. (5.4.6) and (5.4.20)) apply in general, that is, for side charging. When 
constrained to terminal charging the conditions of equations (4. 1.26) and (4. 1 .27) apply. Then the 
linearity requirement for the fractional integral becomes 

- T)'-' (*/{t)+ * j(' - ( f ■ <5A21 > 

and linearity is seen to hold for terminal charging. For the fractional derivative under terminal 
charging conditions, the results of equation (4.1.27) are applied to both sides of equation (5.4.20) 
to yield 0=0, and again linearity holds for terminal charging. 

5.5 Criterion 5: Composition 

Composition or the index law requires proving that 

t d;\ d / /(?)= ,d; , d; /(/)= ,.d; +v /(r) , t>c. (5.5.1) 

The following cases apply: 

Part 1 : the fractional integral case u and v < 0 , 

Part 2: the tools needed for remaining parts, specifically 

a D?f{t)= a D', a Dr'f(t)= a Dr ] „ £>,7(0 for 1 > <? > 0 and 

D"‘ D~ p f{t) = =, D~ P D"'f{t) for p > 0 and m > 1 and an integer, 

Part 3: the fractional derivative case u and v > 0 , and 

Part 4: the mixed case with u < 0 and v > 0. These will be discussed below. The 
integral case will be considered first. 
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5 . 5.7 Part 1, Fractional Integration 


For this part it is required to show that 

c d;\ D;' f(t) = ( D; H -' f(t)= for «>o,v>(U >c. (5.5.1.1) 

The left-hand side then can be written as 

, D ;\.D;'f(t) = j (/ - t)" -1 yy^yj (r - r, )"' /( t, )dr t dr + 

-py-j J {t ~ t)“" ' ift{f -v,a,c,T)dT + yt(h ,-u,a,c,t), (5.5.1.2) 

where h(t) is the input to (argument of) d~ u that would be active during the time period a <t <c. 
Namely, since the i/r's are defined to be zero fora <t<c, h(t ) = a d~'f(t ) for a <t <c (refer 
to figure 5-5-1). 

The case of terminal charging will be considered first; here, for t >c, 

j c 

V r(f-v,a,c,t) = yyyy \{t ~ r)" f{r)dr 

and 

1 c j f 

v(h-u,a,c, t ) = — J(f - r)"' 1 -pypy J (t - T, )"' /( T, )dt l dt. (5.5.13) 


i//(/,-v,a,c,r) \ff(h,-u,a,c,t) 


fit). 



c D;rD;'f(t) 

t > c 


(a) During normal function 



(b) During terminal charging - initialization 


Figure 5-5-1. — Mathematical block diagram composition of fractional integrals (eq. (5.5. 1.1) 
and (5. 5. 1.2)). (a) During normal function, (b) During terminal charging-initialization. 
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Now, for /(/ ) continuous, the first integral in equation (5.5. 1.2) is the uninitialized integral 
product. This is shown to be equal to the uninitialized composition by application of Dirichlet’s 
formula (see Ross (1974); also repeated in Appendix A). Thus equation (5.5. 1.2) becomes 

, d;“ , d;' f(t)= t dr' f(t)+ - ( -m f j (' - 7 )""' ( T - )"' /(*. ) dT < dr + 

W);; 

f f (> - t)""' (t - T, ) v_l /( t, Wt, dr . (55.14) 

r(«)r(v>)|J/ } 1 l) jyi) v ; 

Now working with the RHS and momentarily dropping the common integrands and 1 / TT ’s 

= r d; u ~ v f(t) + J J( )dr l dr + j }( )dz ] dz. (5.5.1.5) 

c a a a 

As long as the integrand is unchanged over the appropriate ranges, the first integral may be 
partitioned as 


= A“~7(') + jj( )dr l dT+]\( )dz x dz+]\{ )dz x dz. 

c a it a u 

Now the first and third integrals can be combined to give 

= t d; u " f(t) + J J( )dz x dz+\]( )dz x dz. 


(5.5.1.6) 


(5.5.1. 7) 


Replacing the integrands, the 1 / TT ’s, and switching limits gives 

,d;\ d;' /(<)=, d - / -t, )' 1 /(r, )dz,dz- 

— j J (r - - r, )' ■' /(r, )dr,dt. (55.1.8) 

r(«)r(v ) J \ 

Again, for f(t) continuous, the Dirichlet formula and transformation may be applied as done 
above (and in Appendix A). Thus, the double integrals become 

, or. DTf(t) = j(' - * /( r -pTvi^ “ r *"*" /(r )l,r 


=,<C"'/M+ 


r(«+v) 


r(«+v)J 

J(t-T)“ + ' '/( T ) dT = t d;“~ r f(t)+i{/(f-u-v,a.c,t) 

a 

r D;\.D;'f(t)= D;"-'f(t), t>c, (55.1.9) 


thus, completing the proof for composition of initialized fractional integrals with terminal charging. 
When side charging conditions apply, if composition is required, that is, 

c d;\.d;' f(t ) , t>c, (5.5.1.10) 
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then 


, d~, u V f(t)+V(f -u-v,a,c,t) = 


1 r, 


I i t 

J( / - t) ^ J( r - rj r J Jt, + — J( ? - t) f-v, a.c.x)dT+ y^k-u, a, c, /), (55L11) 


m- 




rt«) 


for t > c. Applying the Dirichlet formula and transformation (for continuous / (t )) shows the 
leading terms of both sides to be equal. This yields the following relation, which is required to satisfy 
equation (5.5.1.10) 


2 * 

i// (f-u-v,a,c,t) = Y^ij](t-^y , w(f ~v,a,c,t)dT+y/(h -u,a,c,t) 't>c. (5.5.1.12) 


5.5.2 Part 2, Mixed Generalized Integration and Integer Order Differentiation 


The objective of this part is to prove that 

, £>;■, D, ' 7(/)= D, ■/((), 


t>c. 


(5.5.2.1) 


for p > 0 and for m an integer m > 1. This is required to establish composition for the derivative. 

Prior to the main goal of this part the following preliminary result required below is 
established. Consider the following derivative, for 1 > q > 0 , and / > a. 


, D?f{t)= a D) a D^fit) =- D^f{t ) , 

at 


by definition of the derivative. Now 




(5.5.2.2) 


(5.5.2.3) 


by the definition of the fractional integral. Assuming f{t ) differentiable and applying integration 
by parts yields 


( A7(f) s 


d_ 

dt 


+ __L_J ( t _ T ) <-* f '( T ) dr 

dt l r(2-q) T{2-q) { K 

ifzapjia) 

T{2-q) 


■+ u D; {2 -‘ ,} f'(t) 


= D\D; [1 - q] D'f{t)+ 


By part 1, equation (5.5. 1.9), 

„D,'/(r) = 0 , „D 1 - | „D,- M D7(() + 


(I 

T(2-q) 

(t-°V 7(«) 


(5.5.2 .4) 

(5.5.2.5) 
(5.52.6) 


r{i-q) 


(55.2.1) 
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(5.5.18) 


Then, for 1 >q > 0 and for f(t ) differentiable, equation (5.5.2.2) may be written as 

. AVW . K' D)f(t) + ( '~ r ^_^ (n) 


Under the condition that f(t) = 0 V t < a the last term is zero and for 1 > q > 0, t > a 


, d;'/M=„d! . DT'f (<)= .or' . A'/W 


(552.9) 


as the final result. Note, this equation holds from left to right, but right to left fails when 
D'f(t) = 0 or f{t) = constant. This completes the preliminary aspect. 

Again, to establish composition for the derivative it is required to prove equation (5.5.2. 1 ), 
that is, the mixed case of integer order differentiation and fractional integration 

,D; c D; p f{t)= t D; p i D' l "f{t), t>c, (55.U0) 


for p > Oand for m an integer m> 1. 

This requires that 

c d"( c d; p f(t)+\j/{f-p,a,c,t))+y/{h,m,a,c,t)= 

( d; p ( c d?f(t)+y/(f,m,a,c,t))+yf(g,-p,a,c,t), for t>c. 


(5.5.2.1 1) 


where h{t)= tl D; p f(t) and g(r)= a ) . For the case of terminal charging, it has been shown 
that y/(h, m,a, c, t) = 0 and y/(f, m,a,c , r) = 0 for the integer differentiations. Substituting the 
terminal charging ty/ ’s for the fractional integrals gives 

J(f - f(T)dT+j^jj(t - T f(t)dr 


df 


r W 


t(p) j 


1 i /7 m 

for r>c. (55.2.12) 


Now, combining integrals of like integrands 

im f < I ^ 


df 


, t \ t t tm 

(/ - t)'" 1 f( T ) dz j = -^p)l <' " r), '‘ 7r f ^ )dr ' for,>c <5 ' 5 ' 2I3) 


or 

,d; \D;" m= „D7,D';m. < 5 - 5 - 2l4 » 

This conforms precisely to the uninitialized case, which now must be proved to establish the 
initialized case. 

The proof is started with the LHS of equation (5.5.2. 14). Since for the uninitialized case all 
yr's= Othis may be written as 

7(rVf. (5.5.2.15) 
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Integrating by parts m times gives 


= ,D" 


h r( P +*+i) + - D ' f w 


By definition /(/) = 0 V t<a thus f (k) (a)= 0 for all k, therefore 

=„d: M 


(5.5.2.16) 


(5.5.2.17) 


since p > 0, m > 1 , then p + m > m > 1 . Let g{t)= a D"‘ f(t) be a non-zero continuous function of t . 
Then since the integral limits a and t are differentiable functions of t, and since t - 1 > 0 and 
p + m - 1 > 0 then ( t - r)'' +m ’' g(r) and its partial derivative with respect to t are continuous in 
both t and r- Thus Leibnitz’s rule may be applied to give 


,d;\d, 


-(/>+«>) 


# m-\ 


g(t)= 


dt K 


T(^m) T g{t)iT + {l ~ xr "" * (t L 


d "" ' 

dr -■ fl 


D; (p+m -' ] g{t). 


(5.5.2.18) 


Under similar conditions Leibnitz’s rule may be applied a total of m times. The non-integral term 
will always be 0 by the above argument, thus 

. K’M D, -->*(/)= Dr „ />,-/(<) (55.2.1 9) 


d; 


dt 


proving equations (5.5.2.14) and thus, equation (5.5.2. 1) for t >c,by definition of c D q t f(t). 

To establish composition of mixed differintegrals it is also required to prove that for t > c 

, D "\ D; r f{t) = c DT P f{t) (5.5.2.20a) and c D; p c D; /(r)=, D^” f{t) . {55.220b) 

For equation (5.5.2.20a) two cases are considered: 


Case 1 : p> m 

Then under the conditions of terminal charging the LHS of equation (5.5.2.20a) becomes 


d: 


d;'7(0=— 

' J w dt" 


L r (p)f 


j (I - T)'" (' - *>"" 


(5.5.2.21) 


Now, p > m > 1, therefore, p - 1 > 0 and since t - 1 >0, the integrand ( t - r)''" 1 /( r ) and its 
partial derivative with respect to t is continuous in both t and t, provided f(t ) is continuous. 
Then Leibnitz’s Rule may be applied to the combined integral of equation (5.5.2.21) to give 

= ^777 T )"' f( T ) dr (5.5.2.22) 

a 


or 


d'"- ] (p- 1) 

dr~' r(p) 


J(t-Tr ! /(T)dr 


d m -‘ 

dt m -' 


D; {p ~ t] f(t). 


(5.5.2.23) 
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Since p > m, continuity will be satisfied and this process may be repeated for a total of m times, yielding 

=-jT. Dr l ' , — ) /(0=„D, (5.5.2.24) 
at 

Now for / > c , fl D"'~ n f{t)= c D”'- p f(t ) , which completes the case 1 ( p > m) part of the proof. 

Case 2: p < m 
For terminal charging 

j n\ 1 f 1 V 

.. d; , D ;"f( J (* ~ *r l fW'+Tffifc ~ T > f{T)dT \ • (5 - 5 - Z25) 

Now assume m - p > 1 , (i.e., derivative of higher order than integral). Further, assume there 
exists some generalized derivative of order v such that v = n-q = m-p , namely, 

f D;f(t) = d; , d;v( 0. ? > c, ( 5 . 5 . 2 . 26 ) 

where n is the least integer greater than ^ and 0 < q < 1. Now, under the same continuity 
requirements as above Leibnitz’s rule may be applied to the RHS of equation (5.5.2.25) to obtain 

, D: ( /(tVt (5.5.2.27) 

j m - 1 

= — r D Ap -' ] f(t). (5.5.2.28) 

dt m 

Repeat the process b total times, where 0 < p - b < 1, to obtain 

i m-h 

= — — r- a D~ i ' p ~ b) f(t). (5.5.2.29) 

at J \ / 

Let q- p — b then m — b = m + q - p = n yielding 

= a D; u D;‘ l f(t)= i p;f(t) = i D;f(t), frrt>c, ( 5 . 5 . 2 . 30 ) 

thus proving the existence of the derivative and the m — p > 1 case. 

Now, if tn — p < 1, then m — 1 < /?, but by hypothesis for case 2, 0 < p < m thus 

m-\< p<m , 

hence m is the least integer greater than p therefore , D~ p f{ t) is by definition the 

generalized derivative c DJf(t) . Thus, 

l .D' l " c D; p f{t) = c Dr P f(t) (5.5.2.31) 

is proven, completing consideration of equation (5.5.2.20a). 

Equation (5.5.2.20b) is now considered 

, d; p c d;y(/)= <5.5.2.32) 
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Then for terminal charging and t >c> a. 


, D/, D,'7(<K D;\D?f(,) = -£r/(r)rfT. (55.233) 


dr 1 


As above, two cases are considered. 


Case 1: p > m > 1 

Applying integration by parts, for f(t) being m-times differentiable and / <ml * 0 , gives 


(t-rr ] d-'rtr) 


r (p) dr 1 "-' 

Since p > 1, this evaluates to 




r (p) l 




/flj-l 


dr 


m - 1 


/(r)Jr. 


(5.5.2.34) 


('-«r 


r(p) 




-(/>- 1). 


/ m-l 


dt" 


fit). 


(5.5.2.35) 


For /(/), m— times differentiable, and since p > m > 1, the process may be successfully repeated 
a total of m times to get 

f 

k = \ 




r {p-k-i) 

Since by definition /(/) = 0 V t<a, thus f {m ~ k) {a)= 0 for all k, and case 1 (/?>m>l)is 
complete. 


(5.5.2.36) 


fiv)dT. (5.5.2.37) 


Case 2: m > p> 0 and m > 1 

Again, for terminal charging and t > c> a, 

| t j m 

, D;’ . Drf(i)=„D;\ D°'f(t) = 

Again, for f{t) being m-times differentiable and / (m) ^ 0 , integration by parts may be 
successfully applied k times where k is an integer such that k + 1 > p> k> 1, that is, p - k > 0 
then. 


, D;\D'f(i) r( P -j- i) <a> 


(5.5.2.38) 


and because /(/) = 0, V t < a, the summation terms of equation (5.5.2.38) analysis are zero. 

Now 1 > p - k > 0, and since m> p, m - k > 0 and is an integer. Let r = p-k, then the integer 
order derivative can be rewritten as 

.D^'Orfith'Dr'Di'Dr'-'/it). ( 5 . 5 . 2 . 39 ) 
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Then, since 1 > r > 0, using the preliminary result, equation (5 .5.2.9) gives, provided ( , D"' k ' f(t)* 0, 


. D;\ D:'f{t)= u D\D; r u Dr k -'f(t) (5-5.2.40) 

= D 1 D r D l D"'~ k ~ 2 fit) integer order derivatives (5.5.2.41) 

a t a t a t a t J \ / 

=„£>,' ,K .Dr*' 2 /*') by<55.2.9) (5.5.2.42) 

= D] a D~ r a D"‘- k ~ 2 f{t) integer order derivatives. (5.5.2.43) 

Repeating the last four steps a total of m-k times gives 

_ D m-k D~rf( t ) (5.5.2.44) 

il t il t J v 7 

= a D"'~ k ~ r f(t) by (5.5.23 1) (5.5.2.45) 

= u D"'~ p f{t ) , t>c. (5.5.2.46) 


Since u D™~ p f{t)= c D"~ p f(t), for t > c , then equation (5.5.2.32) is true and equation (5.5.2.20b) 
is validated. This completes the terminal charging case. 

For the more general case of side charging the equation 
f d?( c d; p f(t) + y(f,-p,a,c,t))+yf(h,m,a,c,t)= 

r d~ p ( e d"f(t) + Y{f,m,a,c,t)) + v{g-p,a,c,t) for t>c (5.5.2.47) 


must be used directly as the authors have not found a stronger condition as of this writing. 

Some special cases may be considered here. Case A is terminal charging of the fractional 
integrators with side charging of the differentiators. After appropriate substitutions and 
combining integrals in equation (5.5.2.47) 

<T_ 
dt n 


-—r-T J (t - r)'’"' f(r)dT + y/(h, m,a,c,t)= 
r {P)t 

j(f - z) p '' f {m] (r)dr + J(f - t)"' 1 y/{f ,m,a,c,z)dz , 


t > c. 


(5.5.2.48) 


Then with the assumption that f {k] (a) = 0 for all k, the following requirement is obtained 
y/(h, m,a,c,t)= a d; p c d”f{t)- r d™ a d^ p f{t)+ a d~ p y/(f ,m,a,c,t) = 

d~ p \f/(f,m,a,c,t), t>c. (5.5.2.49) 


For the reverse case B, which is terminal charging of the differentiators and side charging of the 
fractional integrators, the requirement is 

j m J m 

\l/(g-p,a,cj)=—y/{f-p,a,c,t)+— ,.d; p f(t)~ / (m) (0 = 

/J " 1 

t > c. (55.250) 
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5.5.3 Part 3, Fractional Differentiation 


The requirement for fractional differentiation is to show, for t > c, 

, D: f D; fit) = D] c D‘;f{t) = f Dr fit), for u > 0, v > 0. (5533) 

Then by definition of the derivative and setting u=m - p and v = n-q with m and n the least 
integers greater than p and q respectively and with 0 < p< 1, and 0 < q < 1 we have 

D;, D'f(t) = Dr At). ( 5 . 53 . 2 ) 

By definition of the derivative 

, A" , d; f{t) = , d; p r d; r d,-*/(0 (5533) 

= D’;\ g (t) . (5.5.34) 

where g(t) = ( .D” V /(/). From part 2 for terminal charging and since g(f) = 0, V t < a, the 
— p andw differintegrals may be interchanged, thus 

< d? , v; f{t) = d; 1 c d; , D;" g (t) (5 535) 

= c D™ +n c D; [p+q) f{t). (5. 5.3.6) 

If p + q< 1, then the definition of the derivative is satisfied and 

c D! c D;f(t)=.Drf(t). (5.5.3. 7) 

However, by definition of the derivatives in equation (5.5.2.16), p + q< 2, therefore if 
p + q> 1, then for f(t) continuous, (/ — t)' 1 /( T ) and its partial derivative with respect to 
t are continuous and we may apply Leibnitz’s rule to the integral to give 

, Dr-' , D; (p * q -' ] f{t)= c Dr fit) , / > C , (5.53.8) 

since the definition of the derivative now applies. Thus, for the terminal charging condition, for 
differintegrals initialized at t = c, the index law holds true for both differentiation and 
integration. 

For the more general case of side charging, the most general situation allows all of the y/ 's 
to be arbitrary. Therefore to satisfy 

,A\-A7(0 = r A" +v /('). for m > 0, v > 0 (5.53.9) 

with D“ +v /(r)= , d‘‘ +v f(t)+y/{f,u + v) requires the following relationship: 
d" +v f{t)+\f/ (/,«+ v)= 

d m tr p d m d-*f(t) + d m d- p d'yf(f,-q)+ 

d"‘d~ p y/(h,n)+ d'"y/{j,-p)+y/(l,m) for t > c, (5.53.10) 
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(5.53.1 1) 


where all i/z's are arbitrary and/i(t), j{t), and l{t ) are defined as 

h(t)= c D; q f(t), 

j(tb r D; c D?M 

and l(th e D;' e D; r D?M 

This is shown in block diagram form in figure 5-5-2. From the results of equations (5.53.7) and 
(5.53.8), specialized to the uninitialized case the first terms on each side of equation (5.5.3.10) 
are equal thus the general side charging requirement becomes 

yf(f,u+v) = d m d- p dy(fr-q)+d m d-y(h,n)+d m yf(j,-p)+yf(l,rn) for t>c (5.5.3.12) 

where the c,t subscripts and a,c,t parts of the arguments have been dropped for convenience. 

Again, special cases may be considered. Terminal charging of the fractional integrators and 
side charging of the differentiators, allows y/(f-q) = a d~ q f(t) and 

yr(j,—p)= a d~ p { a d~ q f(t)+\if(h,n)^to be substituted into equation (5.5.3.12). Similarly, terminal 
charging of the differentiators and side charging of the fractional integrators allows 
yz(/,m,a,c,r) = 0 and y/(h,n,a,c,t) = 0,in equations (5.5.3.10) and (5.5.3.12) above. 



D 

(b) Equivalent diagram 


Figure 5-5-2. — Mathematical block diagram for composition of fractional derivatives 
(eq. (5.53.9) and (5.5.3.10)). (a) Detailed diagram, (b) Equivalent diagram. 
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5.5.4 Part 4, Mixed Integration and Differentiation 


It still remains to prove the general mixed case 

, D * , D;'f(t)=D;' c D';f(t)=D'rf(t), t> C , (5.5.4.1) 

for u > 0 and v > 0. 

For the case of terminal charging, let u = m— p where 0 < p < 1 and m> 1 and an integer, then 
using part 1 results the LHS becomes 

r A" c D; v f(tb c D? c d ; 1 ’ c D; r f(t)= c d; t . z>r (,,+l) / (')• (5.54.2) 

Now, under the conditions allowing part 2, this becomes 

= .D; {p+r \.D;f(t) (55.43) 

— C D^' c D~ p c D™ f (r) (5.5.44) 

= c D;\,D"f{t). (55.45) 

This proves the first equality in equation (5.54.1). Now, from equation (5.54.3) with v + p > 0 
therefore, using part 2 results this gives, 

, A", D; v f(t)= c Dr iv+p) f{th c Drf(t), t > c , (5.54.6) 

completing part 4. This completes criterion 5. 

5.5.5 Comments on Criteria 

The totality of parts 1 , 3 and 4 of criterion 5 above combine to give the following result. 
Under the conditions of terminal charging of the u th and v th differintegrals 

, d : , d ;/(0 = d;\ d;/(0 = ,zr v /(0 for r > c (5.5.5.1) 

under the following conditions: 

(1) « < 0 , v<0 for f(t) continuous 

(2) w > 0 , v > 0 f(t) is m -times differentiable, a D'"f(t) exists and is 

a non-zero continuous function of t, t > a ,where m an 
integer is the larger of lnt[u] or Int[v] . 

(3) u < 0 , v > 0 same as (2). 

We now discuss the basic criteria results derived above for a generalized (fractional) 
calculus. The scope of this work was limited to the real domain. Attention will first be focused 
on the terminal charging case. 
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Terminal Charging 

Relative to criterion 2, backward compatibility with the integer order calculus, the addition of 
the initialization function is clearly a generalization relative to the integer order calculus. In a 
strict sense t j/(t) * 0 violates criterion 2, however we are seeking a generalization of the integer 
order calculus and it is clear that this generalization (i.e., the formal addition of an initialization 
function) will be very useful in many applications. 

Relative to criterion 3 the zero order property holds for terminal charging. 

Relative to criterion 4 linearity holds for the terminally charged case. 

Relative to criterion 5, composition holds for terminal charging subject to the constraints 
noted in equation (5.5.5. 1). It is noted that the constraint, / (t) (c) = 0 Vfc,no longer exists, which 
satisfies a primary objective of this work. This constraint has effectively been contained in 
(shifted to) the requirement /(/) = 0, V t < a . This allows practical initialization for fractional 
differential equations. 

In summary , the terminal charging case is backward compatible with the integer order 
calculus and satisfies the applicable criteria established by Ross for a fractional calculus. 

Side Charging 

The state of affairs with regard to side charging is less definitive. Criterion 2, backward 
compatibility, is the same as the terminal charging case since the i/r ’s are arbitrary. 

Relative to the zero property criterion 3, the conditions that 

\jf(f -p,a,c,t) = ^—\{t-TY and y/(h,m,a,c,t) = 0, (55.52) 

are not required for side charging. When these conditions are not met, the zero order operation on 
/(f) will return f(t) + g(f) , that is, the original function with an “extra time function ( gU)), 
the effect of the initializations. At first reading this may appear somewhat bizarre, however, this 
is of benefit (realizable) in analyzing certain physical systems. This will be demonstrated in the 
applications section of this paper (see for example the generalized zero order operator), and the 
reader is requested to reserve judgment until that point. 

Relative to criterion 4, linearity for the side charging case, demands the requirements of 
equation (5.4.6) or equations (5.4.18) and/or (5.4.20). 

Relative to criterion 5, composition for the side charged case additionally requires the 
satisfaction of equation (5.5. 1.1. 2) for fractional integration, equation (5.5.2.47) for mixed 
generalized integration and integer order differentiation, and equation (5.5.3. 12) for generalized 
differentiation. It is noted that with equality of the leading terms on the two sides of equation 
(5.5.2.47) the requirement is simplified to 

, </;>(/ -p)+v{h,m)= r d; p ii/(f,m)+\]/(g-p). (5.5.5.3) 
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These are not so much of an issue as first appears for practical application. In the solution, for 
example, of fractional differential equations, y/ (t) will be chosen in much the same manner as 
initial conditions are currently chosen for ODE’s or for PDE’s. This will imply the nature of 
fit) fora <t<c. The new aspect is that to achieve a particular initialization for a given 
composition now requires attention to the initializations of the composing elements. That is, the 
requirements derived above must be used to determine the new initialization. 


6. Laplace Transforms of Generalized Differintegrals 

In this section general forms will be derived for the Laplace transform of fractional 
derivatives and integrals. The derived forms should be generalizations of those available for use 
with ordinary (integer) differential equations and will include the effects of the initialization 
functions for the fractional differintegrals. 


6.1 Preliminaries 

It is useful to consider the following two simple cases before attempting general results. 

Based on the integer order calculus, when /(f) is piecewise regular and of exponential order, the 

t 

Laplace transform of J f{t)dt is given by 

b 

l\ j = - L{f(t)} + - j f(t)dt (6.1.1 ) 


(see e.g., Wylie (1975) p. 265). Contrast this to the generalized calculus case for the same 
integral. 


L{ * D;' f(t)} = L j f(t)dt + y/(f-\,a,b,t)\. 


( 6 . 1 . 2 ) 


If equation (6. 1 . 1 ) is used to evaluate the transform, we have 


1 1 ^ 

4. A '/(')} = - 4/(0} + - J 10)41 + 4 V (/,-!, a, b, 0 }. (6.13) 


However, in the most general case y/ is arbitrary, if it is chosen as yf = const then since 

, r , const 
L{const\ = , 


(6.1.4) 
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it is clear that this term contains the initializing effect of the second term on the RHS of equation 
(6.1.3). Hence it is not necessary to include such terms that redundantly bring in the effect of 
initialization from the integer calculus. Thus, 

L{ h D,~‘ f{t)} = ^ L{f(t)} + L{w (/ -1, a, b, /)}. ( 6.1.5) 

Consider now the integer derivative. Based on the integer order calculus under appropriate 
conditions (see Wylie (1975) p. 264) the Laplace transform of /'(/) is given by 

L{f'(t)} = sL{f(tj}-f( 0 + ). (6.1.6) 

Contrast this to the generalized calculus case for the same derivative 

t{«D, '/(<)} = i{^ + V'(/,l.«,0,r)J • (6J.7) 

If equation (6.1.6) is used to evaluate the derivative in equation (6.1.7), the following is obtained: 

L{ 0 D'j(t)} = sL{f(t)}-f(0 + ) + L{y/(fXa,0,t)}. (6.1.8) 

Again in the most general case y is arbitrary, if it is chosen as ip = -f(t)S(t -0 + ), then, since 

l{-/W5(/-0 + )} = -/(° + ) (6.1.9) 

(see e.g.. Gabel and Roberts (1973) pp. 72-77) it is clear that this term contains the initializing 
effect brought in by the integer order calculus term. Again it is not necessary to include such 
redundant terms. 

The redundant terms are introduced into the above equations by the product term of the 
integration by parts formula, which is applied in deriving the reference equations in the integer 
order calculus. In the subsections that follow, these redundancies will be eliminated without 
hesitation when it is appropriate to do so. 

In the material that follows, considering discontinuities at t = 0 , evaluations stated as t - 0 
will be understood to be evaluated at t = 0 + , consistent with the conventional definition for the 
Laplace transform (Wylie (1975)). 

6.2 Laplace Transform of Fractional Integrals 

For simplicity the starting point is taken as c = 0 . Thus it is desired to evaluate 

l{„D, "«/(')} - f{r)dr + \i/{f-q,afi,t) dt q> 0, r>0. (6.2.1) 
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The convolution theorem for the Laplace transform (Churchill (1958)) may be applied if 
H(s ) and G(s) are transforms of h{t) and g(t), which are sectionally continuous functions and are 

of exponential order as r — » °o . Then the convolution h(t)* g(t ) exists and is defined as 

1 

h{t)* g{t) = j h(r)g(t-z)dT (6.2.2) 

o 


and 


L{h(t)* g(t)) = H(s)G(s) = /I J h(r)g(t - t)dr 


To 


(6.2.3) 


t<i- 1 


Then taking h(t) = f{t) and g(t)= — tt the convolution theorem (for eq. (6.2. 1 ))yields 

rw) 


( 1 A 


F(s)G(s) = L(f{t))L r = — L(f(t)) q> 0. 

l r W)J 


(6.2.4) 


This yields the general result, 


(6.2.5) 


Now in equation (6.2.5), t/r(/ ,-q,a,0,t) may be thought of as the composed (equivalent) 
effect of the initializations of the mathematical elements used to create 0 D~‘‘ f(t) . For the 
fractional order case there are infinitely many ways 0 D^f(t) may be composed, as opposed to 
the integer order calculus case where only combinations of integer order integrations are 
possible. For example, in the integer order calculus, the familiar Laplace transform. 


L{f in) (t)} = s n F(s)-s"-'f(0+)- S n - 2 f'(0+).... 0+), (6.2.6) 


infers that /* " 1 (t) is composed of (or decomposed into) n separate differentiations each of order 1 . 


6.2.1 Integer Order Decomposition of Fractional Integral 

It is clear that the composition law must be satisfied to achieve equivalent initialization 
functions. Further, so long as only generalized integration is considered, composition is satisfied 
(eq. (5.5.5. 1)). 

As an example decomposition, the following is obtained through integration by parts; for 
q > 1, consider 

^M/ a A 0) = j e ~' ! w(f ~q,a, 0, t) dt . ( 6.2.1.1 ) 

0 
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Now integrating by parts, with 


u = if/(f-q,a,0,t). 


dv = e- s 'dt. 


then 


and with —j— = I//” yields 
dt 


du = — lj / (f -q, a ,0, t)dt , 


v = - 


— S 


( e” 


L(y/(f -q,a,0,t)) = —y/{f-q,a,Od) 


( 6 . 2 . 1 . 2 ) 


( 6 . 2 . 13 ) 


+ - \e x '\f/ {U (f-q,a,0,t)dt (6.2.1.4) 

f J 


or 


L(y/(f ,-q,a,0,t))=0 + -y/(f ,-q,a,0,t\ =o + ~L(if/ { ' ] (f ,-q,a,0,t)). (6.2.1.5) 


Repeating this process gives 

l[\f{f,-q,a&t)) = ^ + “( V'(.A^«A')L + ^ j. ( 621 . 6 ) 

Now repeating the process a total of n times, where n is an integer such that n + 1 > q > n gives 

/.(^(/,- 9 , O ,O.0)=^i)r , ’ , (/.-«.‘>.O./)) + S^V' IJ - ,, (/.-?.a.O.<)|, - , (6.21.7) 

j = 1 A 


as the expression for the equivalent initialization function. This equation then yields 



( 6 . 2 . 1 . 8 ) 


The inference of this equation is that the q th differintegral is composed of n order 1 integer 
integrations and a fractional integration of order q-n (see fig. 6-2-1-1). Further, the order 1 

integrations are each initialized by a constant (if/*' 1 /=0 terms in the summation). This form is 

useful in showing backward compatibility with contemporary Laplace transform theorems for 
repeated integer integrals in the integer order calculus, namely for n integrations (see Appendix B), 


j / >i t „_ i 

L |JJI -• 

I o o o o 



n = 1,2,3.... (6.2.1.9) 


The compatibility is seen by taking q = n = 1,2,3... and properly selecting if/ in equation 
(6.2. 1.8). In equation (6.2. 1.8) the summation term is the result of repeated application of 
integration by parts and the term represents a redundancy relative to the preceding term for the 
generalized calculus. 
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(b) Laplace transform n -° rder 1 Integrations 


Figure 6-2- 1 - 1 . — Block diagrams for integer order decomposition of the fractional integral (eq. (6.2. 1 .8)). 
(a) Mathematical diagram - time domain, (b) Laplace transform - frequency domain. 

It is clear that infinitely many more possibilities exist than this form. As a further example the 
integer order integrations in figure 6-2- 1-1 could just as well be replaced by generalized 
(fractional ) integrations each also of order 1. Each of these would allow a non-constant 
initialization function; thus consider the decomposition shown in figure 6-2- 1-2. For this case, 
o D~ q f{t) is decomposed into n generalized order 1 integrations and a fractional integration of 
order - q + n . Thus, n is the greatest integer less than q. Since all of the operations are 
integration, composition holds, and the following may be written relative to this diagram: 

oA"V(0=V^+o^rVi(0 . q>0, 

= ¥x +od;'(yr 2 + 0 4"V 2 (0) . 

= ¥i + o d; l yr 2 + 0 d; 2 (y/ i + 0 d;'x n _,(t ) ) , 

= ¥,+ 0 <CV 2 + o d; 2 ¥i +o 4" V 3 (0- (6.2.1.10) 


y/(f,n-q,a,0,t) r) V',(jr„. l .-l.a.O,/) 

I — 1.4 . r~ 1 4 r— 14 j — hi i — 1 4 

I L_ I + v : L_ I '• L_ I + x ’> 1 L 1 I + x ^ D < 


n-Generalized Order 1 Integrations 


Figure 6-2- 1-2. — Mathematical diagram for generalized integer order decomposition of the 
fractional integral of order q (eq. (6.2.2. 1 2)) for t > c = 0. 
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Repeating this process a total of n times and observing that 

x, 0) = + n,a, 0, t ) + 0 d;‘ l+ "f(t ) , 


( 6 . 2 . 1 . 1 1 ) 


gives the following: 

o d;v( 0= o d; q f(t)+ o4~>(/ ,-9+«,a,o,/)+Io^ (r 'V;h ,-u,o,f), (6.2.1.12) 

H 

as the mathematical representation of the diagram. The Laplace transform is now taken of this 
equation while noting that the effects of initialization are explicitly included and should not be 
included again when transforming the uninitialized integral terms. 

Thus, 

(621,3) 

Because L{const } = const / 5 , it is clear that the difference between this equation and equation 

(6.2. 1.8) is the effect of initializing here by a function of time instead of a constant as in equation 

(6.2. 1.8) . This equation (6.2.1.13), of course, is a generalization of equation (6.2. 1.8). Further, if 
q is an integer, namely g = n,then xft) = /(?),and \//{f ,—q + n,a,0,t) may not be required. 

Then this equation specializes back to the form for the Laplace transform of a multiple integer 
integral from contemporary transform theory (see Appendix B), namely, 

» = 1.2,3... (6.2.1.14) 


6.2.2 Fractional Order Decomposition of Fractional Integral 

A much more general decomposition of 0 D; q f(t) can be obtained that is not limited to integer 
integral elements. Consider the mathematics associated with the decomposition indicated in figure 
6-2-2- 1 . Since only fractional integrations are considered, composition holds for the case shown in 
this diagram. Thus, dropping the subscripts on 0 for convenience we have 

jc,(/) = jc„ +l (f)= a D; q f(t) , t> 0, q k > 0 Vfc, (6.2.2.1) 

n 

where q = ^q ( . Then starting from the inside 

i=i 

L{x 2 (/)}= = l{ 0 d; q 1 JC, (0} + l{ iff, (jc, -q , , a,0, t) } , ( 6.22.2 ) 

for which the following shorthand will be used: 

l{jc, } = }-m/z, ( 5 ) . (6.2.23) 
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• • • 


W,0) Vn(‘) 





D h 



~v 

D c{ D^'- 

— v 

n- Fractional Integrations 



Figure 6-2-2- 1 . — Mathematical diagram for multiple generalized (fractional) integrations, t > c 


Then 

and 


l { jc ,} = s ^ + y/ { (s) 

l { jc 3 } = L{d~ q 2 x 2 } + y/ 2 (s) = + \f/ 2 {s) 

= s ~“ 2 { l { x , } + ( 5 )} + y / 2 (s ) ; 


again 

l{x 4 } = L{d^ q 'x 2 ] + i// 3 (s) = + y/ 2 (s) 

= s~ q, s~ q2 s~ q> l{jc,} + s~ q, ~ q '-y/ f {s) + s~ q, y/ 2 (s) + y/ 2 {s) 
= s - q '- q -- q 'L(x ] )+ s^ 2 i//, (.?) + s~ qi y/ 2 (.v) + 1//3 (.v) ■ 

Repeating the process to x n+l 

L{x „ +1 } = s - q '~ q '--- q ° /.{*, } + j-*-*---*- yf l (s) + s - q >- q *-~- q ' y/ 2 (i) + 

• • • + s~ q "'- q « y/ n _ 2 (s) + s- q ° 1 j/ n _ x (5) + y/„ (s) . 

n 

Now defining B a = equation (6.2.2.10) can be written as 

i=a 

^{•*n+l} = lS + 5 B: l//,(5)+ +S Bn il/ n _ x {s) + y/ n (s), 

or more compactly as 

L{ 0 D; q f(t)}=L{ 0 D; B 'x l (t)} = L{x m+ ,}= 5- a ‘L{x,} + + 

7=1 

n 

where B t = q ( . 1 < a < n q k >0 V& . 


,/>r z vto 

= 0 . 

(6.2.2.4) 

(6.2.23) 

(62.2.6) 

( 6 . 22 . 7 ) 

( 6 . 22 . 8 ) 

(6.22.9) 

( 6 . 2 . 2 . 10 ) 

( 6 . 2 . 2 . 11 ) 

(6.2.2.12) 
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The attraction of this form is the fact that the q { can be integer, non-integer or mixed. Further, 
although this analysis was done for q > 0 , that is, for fractional integrals the basic form will also be 
shown to hold for fractional derivatives provided composition is satisfied. The effective initialization 
then for this case is 


l\Weffea } = iM*) + X ^ V J W ‘ 

7=1 


(6.2.2.13) 


6.3 Laplace Transform of Fractional Derivatives 


Again, for simplicity, the starting point is taken as c = 0. Here it is desired to evaluate 

Mo AV(0} = L{ 0 Dr 0 D; p f{t)}, U > 0 , ( 6 . 3 . 1 ) 

where m is the least integer > u such that u = m- p. Then substituting the definitions and 
expanding, gives 


L{ 0 D:f(t)}=L{ 0 D;\D;'’f(t)}=L — [j ^ r ^) f(r)dt 




L|— w(f-p, a, 0,r)| + L{i/r(/i,m,a,0,t)}, t> 0, (6.3.2) 

where h{t)= a d~ p f(t). The sum of the last two terms on the right-hand side can be thought of as 
the equivalent y/ for 0 Z)"/(t) , thus 


d 

y eq (f,u,a,0,t)= — 0, t ) + \jf{h, m, a,0,t) . 


Consider the first term of the right-hand side of equation (6.3.2), 


df 


»>-*)■ 

I r ( P ) 


7-1 


-f(r)dr 


= f e ~ s ' — £—\ 

J dt ' 


d"'~’ 'dt-r) 1 ’- 


dr-'{ r{ P ) 


■f{z)dt 


dt. 


Applying integration by parts where 


( im-l 1 


u=e 


dv=- 


then 


du = -se s ' dt , 


dt 


v- 


=rj 


(t-rT' 


dt m ~ x { r ( P ) 


-f(r)di 


dt 


i m- 1 * 


(t- T)"-‘ 


-f(r)dT , 


yields 


i rti - 1 f 


- e 


(r-r) 

j<-'{ r(p) 


■J 


P-1 


f(r)dT 


df-'{ r( P ) 

( d”'- 2 Ht-r)"- 


+ s I e 


± 

dt 


dt- 2 l r ( P ) 


f(r)dr 


(6.3.3) 


(6.3.4) 


(6.3.5) 


(6.3.6) 


dt. (6.3.7) 
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The first term evaluates to zero (further, this term would be a redundant initialization), thus 


dt' 


f(r-r) 




Vo 


r (p) 


f{t)dT 


= sL 


( jm- 2 I 


dt 


tJ 


(t-tr 


dt m ~ 2 { r(p) 


f(r)dt 


(6.3.8) 


Repeating this process a total of m times yields 

' d m 'c(t-r)"-' 


^l^w nr)dt r LUtfil)l 


(6.3.9) 


Therefore, one form of equation (6.3.2) is; 

l{„d: f(,)} = l{„d: „d;' f(t)} = 


s' L{„d;"f(t)}+ L\ ~ V {) ,-p,a,0,l)\+ I )} (6.3.10) 


or with equation (6.3.3), 

L{ 0 D;j(t)} = s m L{ 0 d; p f(t)} + L{y,{f,u,aAtj\, (6.3.11) 


where i j/ eq (f,u,a, c,t) is written asi f/(f ,u,a,c,t) to allow the generalization that follows. 
Applying the results for the fractional integral, equation (6.2.5), gives 


L{ 0 D‘;f(t)} =s m - p L{f(t)}+L{y/(f,u,a,0,tj\ 
= s u L{f{t)} + L{\ff(f,u,a,0,t)}. 


(6.3.12) 


These simple forms are the most general forms for the Laplace transform of 0 D‘‘f(t) . It is 
noted that this form is the same as that of equation (6.2.5), thus in equation (6.3. 12) u may take 
on any real value. 

Now, as in the fractional integral case, there are an infinite number of ways in which 
0 D,"/(r) may be decomposed. Thus, there are many possible formulations for t ff eq (f,u,a,c,t) . 
Some of these will now be explored. 


6.3.1 Integer Order Decompositions of the Fractional Derivative 


An equivalent form may be derived for L|i/r(/,«,a,0,/)| with a composed of m integer 
order differentiations and a single order p fractional integration. Consider then, the second term 
of the RHS of equation (6.3.2) or (6.3.10), 


j m I 00 i 

L\ — — l if ( f-p , a, 0, r) > = J e st 

J a 


dt r 


dt 


dt' 


-y/(f-p,a,0,t) \dt. (6.3.1.1) 
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Again, integrating by parts with 


im-l 


u - e 


du = —se "dl , 


d d n 

dv = - y(f -p, a frt)dt , 

d w_l 


(6.3.1.2) 

( 6 . 3 . 13 ) 


Therefore, 


' , m t jm-l °° ~ ^">-1 


dt 


lo 0 


/m-1 


dr* 


•VM 


(/,-/?, a,0,t)| +sL 1“pr^(/-P. fl *°*0f- ( 63L4 > 


l/=0 


This result can now be applied to determine 


d d 


m-2 


dt dt 


y/(f,-p,a,0,t)\ = 


dt' 




1=0 


+ sL<-^ry/(f-p,a,Q,t)\. (6.3.1.5) 


Then, 


i/ti-1 


L\ — V ( f~P , a , 0 , t) \ = - y (/ -p,a,0, 0 ] 


dt 


1 = 0 




Jr” 


r v(f,-p*a,0,t) 


t = 0 


d m - 

+ sL\-^ry/{f ,-p,a,0,t)\ 


(6.3.1.6) 


Repeating this process a total of m times and writing 


d k 


dt 


7 V(f-p,a&i) = yr lk) (f-p, a ,0,t), 


(6.3.1.7) 


yields 




7=1 


It is noted here that the summation term, in this equation, is the result of repeated integration 
by parts from the integer calculus. Clearly an arbitrary if/ may be chosen such that the 
initialization information contained in the summation can be incorporated into the previous term. 
That is, the summation term is redundant but will be carried so that an explicit comparison of 
form back to the integer order calculus may be made. 


NASA/TP— 1 998-2084 1 5 


55 



Substituting equation (6.3. 1.8) into equation (6.3.10) gives the following expression for the 
Laplace transform of the generalized derivative (decomposed into integer differentiations with 
V't"' -7 ' (/,-/>, «,c,r)( r=0 = constant initialization): 


L{ i) D:f(t)} = L{ 0 Dr 0 Drf(r)} = 

m 

+ I +L{v(h,m,a,0,t)}. (63.1.9) 

j=] 


| r 

Here, 0 d~ p f (t) = ( t - t )' 1 f(r)dT, is the uninitialized fractional integral originating at 


t = c — Oand h(t)= u d~ p f(t) as defined for equation (6.3.2). Equation (6. 3. 1.9) is a general 
expression for either side charging or terminal charging, (but is limited to the case where the y/ 
terms in the summation evaluate to constants). Equation (6.3. 1.9) infers the mathematical block 
diagram shown in figure 6-3- 1-1. 



Figure 6-3-1- 1. — Mathematical block diagram for equation 6.3.1.10. 


The first term of the RHS may be written as .v"l{/(/)}, and the above equation becomes 


l{ 0 d; /(,)}= l[ 0 d;„d ; '/(<)}= 

m i 

s H a, 0, r)| u + L{y/(h, m,a, 0, / )}. 


(6.3.1.10) 


Note, what has been done here is to decompose the integer order derivative part of 0 D" f(t), 
namely, 0 D"‘ into m integer differentiations, each of order 1 . Specialization of this equation 
back to the integer order calculus expression for the transform of the derivative, namely, 

m 

L{/ w W} = fL{/(f)}-£/ l /^)(0 + ), (6.3.1.11) 

7=1 
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can be seen as follows. In equation (6.3. 1.10), take p = 0 (zero order operator; that this operator 
can have an initialization in the generalized calculus will be demonstrated in the application 
section). Set y/(h,m,a,0,t) = 0 then y/(f,0,a,0,t) can be specialized such that 

= / (m - y) (0 + ) (6.3.1.12) 


and the reversion is shown. 

In the context of equation (6.3. 1. 10), 

L{y/(f,u,a,Oj)} = -^s H (y/ {m ' j) (f-p,a,0,t)^ o + L{y/{h,m,a,0,tj \ . (6.3.1.13) 

j=\ 

Further generalization of equation (6.3. 1 . 10) is possible in several ways. From equation (6.3.2) or 
(6.3.10) consider, again, the term L 


y/(f-p,a,0,t)\. Now since the initialization can be 


completely arbitrary, let 




.. +„<0„M= S .<//■>;(»)■ (63.1.14) 


Then, 



Y J s J L{y/ ] (t)} , 


(6.3.1.15) 


where the redundant terms have been dropped. Then equation (6.3.1.10) can be generalized to 


L{ 0 D;m}=L{„D;,D ; '/(')}= 

m 

s"L{f(t)} + X s i L[\p j (t)} + L{y/{h, m,a, 0, t)}. 


(6.3.1.16) 


The difference in the two equations, of course, is that here the order 1 derivatives are each 
initialized by time varying functions. The inference of this is that the order 1 differentiations are 
now generalized order 1 differentiations. The derivative initialization L^xpih, m, a,0 , /)}can also 
be similarly distributed over the order 1 derivative terms to produce a similar effect. 


6.3.2 Fractional Order Decompositions of the Fractional Derivative 


The Laplace transform of the generalized (fractional) derivative may also be based on 
decomposition into fractional derivatives and integrals. The general fractional order 
decomposition of the fractional integral as performed in a previous subsection may be extended 
to the case with fractional derivatives. To do this, replace the fractional integrations indicated in 
figure 6-2-2- 1 by fractional differentiations; let- q k =r A with r k >0 (refer to fig. 6-3-2- 1). 
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y/(x k -p,a,0,t) y/(h k ,m k ,a,0,t) 



(a) Differentiation element 



(b) Equivalent differentiation element 

Figure 6-3-2- 1 . — Fractional differentiation element, (a) Differentiation element, 
(b) Equivalent differentiation element. 


Then the decomposition is valid provided: (1) that composition is satisfied at each step of the 
reduction and (2) that for each k such that r k > 0, the initialization y/ k is considered as the 
equivalent y/ for differentiation, namely from equation (6.3.3), 


u k 

{* k . r k a,0, t) = —— y/(x k -p k _ x ,a, 0, r) + y/(h k , m k ,a, 0, t) , (63.2.1) 


where m, is the least integer such that r k = m k - p k and h k = tl d l Pl x k {t) . The result then by 
identical mathematics as in the integral case is 

(63.2.2) 


Clearly the extension of this result to all r k additionally requires the full applicability of the 
mixed-mode form of the composition law (eq. (5. 5.4. 8)) and the use of the equivalent y/ as given 
in equation (63.2.1) whenever r k > 0. This general result, equation (63.2.2), for all r k is 
probably the most powerful and useful of the forms derived in this subsection. 



6.4 Effect of Terminal Charging on Laplace Transforms 


In the results generated in this section all y/ terms were treated as completely arbitrary, that 
is, a side charging assumption is implicitly made and the forms above apply to either terminal or 
side charging. The effects of terminal charging may be readily determined by appropriate 
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substitutions for the yr terms. That is, for the integer derivative terms substitute 

\f/(h, m,a,0,t) = 0 m = 1,2, . . . , and for the fractional integrations (either as part of the derivative 

or alone) take 


y/(f,-p,a, 0, t) = ^j } (t - r)" _1 f(r)dr. 


(6.4.1) 


Of course it is not necessary to terminal charge both terms of the derivative, but the choice is 
made based on the desired decomposition. 

Some modest simplifications may be achieved. For example, relative to the Laplace 
transform of the fractional integral, the y/ U ) terms in equation (6.2. 1.8), if terminal charging is 
assumed, are 

M 




t > 0. (6.4.2) 


l ;=0 


Applying Liebnitz’s rule gives 


d‘- ! ( d ?(„-! X>-r) (lHH 




- -f- 

2 Jt J 


dt k dt 


r («) 


f(r)dr 


t> 0 . (6.4.3) 


7=0 


Now because ( q - 1)/ r(y) = 1/ r(q- 1), the integral is recognized as the initialization of the 
(i q - l)th order differintegral and continuing the process yields 

y/ {k) (f-q, a ,Oj)[ =o =y/(f,-q + k,a,0,t\ =o , t> 0. (6.4.4) 

In similar manner, the more general term in equation (6.2. 1.8) yields 

y/ (n) {f-q, a, 0, t) = y/{f-q + «, a,0 , j / ( r ) dr ’ t>0 ' (6A5) 

where n + l> q > n>0 and n is an integer. 

6.5 Effect of Starting Point 


6.5.1 Fractional Integral 

The effect of starting with non-zero a (and c ) is considered for the fractional integral. 
Specifically, the case of c> a> 0, and q > 0 is considered. Then for the uninitialized fractional 
integral 

l {„ d , -*/(')} = i\ a = L {f^) 1 _ r)<H f^ dT \ q>0 - (6 ' 5 ‘ U) 
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(6.5.1.2) 


By definition f(t) = 0, for t < a, therefore we can write 


f{t) = u{t-a)f{t ), 


where u(t) is the unit step function, and thus the result is 

i{- D .‘ , /(')}=/-{„<#;V(')} = 4« rf ." , W'-a)/('))}. y >0, a> 0. (65.13) 


6.5.2 Fractional Derivative 


Under the same conditions as above, the uninitialized fractional derivative is considered, then 

f I m t / ] 

Lid?m}=L{.dr =t|^r-p^J(»-i-r 7 w<(rj , «>o, (6.5.2.D 


and as before m is the least integer > q such that q = m - p . Again by definition, 
f{t) = 0 for t < a , therefore we can write 

f(t) = u(t -a)f(t). 


Thus 


tm 1 1 

c < /(')} = L< — J (/ - a)f(z)dT 


1 


(6.5.2.2) 


(6.5.23) 


= L\— 0 d; p (u(t- a) f(t))\ = a) f(t))}, q> 0. (63.2.4) 


Thus, combining results of equations (6. 5. 1.3) and (6. 5.2.4) gives 

l{a '/(<)} = i{„rf, "(«(>-«)/(»))}. 


(6.5.2.5) 


6.6 Laplace Transform of Initialization Function 


6.6.1 Fractional Integral 


It will sometimes be useful to deal with the initialization function in the Laplace domain. 
From Laplace transform theory (Wylie (1975) p. 281) the so called first shifting theorem is given 
by 

L{g(r-a)w(/-a)} = e" m L{g(f)}, a > 0. (6.611) 
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Letting f(t ) = g(t - a ), this can be rewritten as 

L{f{t)u{t-a)} = e~“'L{f(t + a) }, a> 0. (6.6.1.2) 

Now for the fractional integral (under terminal charging), the desired Laplace transform then is 

L{ V (f.-q, a, c, 0 } = l{ „ d- ■•/(')}- L{ , </,-*/ ( 0 } • (6.6.13) 

Using the above result, equation (6.5. 1.3) 

= L{ 0 J; l 'f{t)u{t-a)}-L{ 0 d; ll f(t)u(t-c)} (6.6. 1.4) 

= s' q L{f{t)u{t - a)} - s- q L{f{t)u{t - c)} (6.6.1.5) 

and applying the modified shifting theorem gives the result, for q > 0, 

,-q,a,c,t)^ = s ~ q \e ~“ s L{f(t + a)}-e ‘ 'L{/(t + c)}], c > a > 0. (6.6.1.6) 


6.6.2 Fractional Derivative 

The equivalent initialization function for the fractional derivative c D“f{t)= c D'" C D, ■ P f(t ) is 
given by (eq. (6.3.3)) 

(1 m 

y/ eq {f,u,a,c,t) = — \l/{f-p,a,c,t) + y{h,m,a,c,t), u> 0, (66.2.1) 

where m is the least integer > u such that u = m - p . For terminal charging it has been show that 
y/(h, m,a, c,t) = 0 , therefore 

L{\i/(f,u,a,c,t)}= s m L{y/{f-p,a,c,tj\, (6.62.2) 

and applying the above result (equation (6.6. 1.6)) for the fractional integral yields 

L{\j/(f,u,a,c,t)}= s u {e'“ s L{f(t + a)}-e' cs L{f{t + c)}}, c>a> 0. (6.6.23) 

6.7 Summary — Comments on Laplace Transforms of Differintegrals 

(1) The general case forms of equations (6.2.5) and (6.3.12) for the fractional integral and 
derivative, respectively, will be most useful in many applications. Subject to satisfaction of 
continuity and composition/index law as required in the derivations, these solutions may be 


combined into the single form 


L{ 0 D?f{t)} = s‘‘L{f{t)} + L{y/(f,q,a,0,t)}, 

\/q. 
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Where for fractional integration (q < 0), \p(f ,q,a,0,t) is arbitrary and for fractional 
differentiation ( q > 0) where q = m - p and m is the least integer > q , then 

cf m 

\j/(f,q,a, 0 , t) = — y/(f -p, a, 0, t) + y/(h, m,a, 0 , t) 

where h{t)= a d~ p f(t). The use of equations (6.7.1) and (6.7.2), along with block diagrams, to 
clearly delineate the component parts of the various elements (the decomposition) can be an 
exceedingly useful approach. 

(2) The solution of integrodifferential equations containing the combination of fractional 
order and integer order terms ( non-generalized ) can effectively apply the forms of equations 
(6.2. 1 .8) and (6.3. 1.10) for initialization by constants for the integer order terms. When the 
forms contain generalized integer order terms the equations (6.2.1.13) and (6.3.1.16) should be 
applied. Of course these later two equations are generalizations of the former two and may with 
care be used in both cases. 


(3) The forms for the Laplace transforms with fractional decompositions are still more 
general than either pair of equations discussed in (2) above. These are equations (6.2.2. 1 2) and 
(6.3.2.2) and are nearly identical in form. In fact these equations taken together may be 
generalized into a single result, which is equation (6. 3.2. 2) with r k free to take on both positive 
and negative values, that is. 


Mo D ’' x i ( r )} = L { x «+ 1 (0} = * B ' M*] (0} + ¥ n (•') + X Vj ( 5 ) ’ 

j = 1 


/I 

where B a = ^ r ; , V r ( , 1 < a < n. 

i=a 


(6.7.2) 


For this equation, when r k >0 it is understood to be interpreted as r k - m k - p k as derived 
in (6.3.2. 2). 


(4) Many of the forms derived in this section on the Laplace transform have been derived to 
demonstrate how the initialization function generalizes the Laplace transforms seen in the integer 
order calculus. An example of this is the reduction of equations (6.3.1.10) and (6.3.1.16) to the 
familiar form of equation (6.3. 1 . 1 1 ) for the nth derivative, when y{t) is specialized to be a 
constant. Similar is the reduction of equations (6.2. 1.8) and (6.2.1.13) to equation (6. 2. 1.9) for 
the Laplace transform of multiple (integer order) integrals. 

(5) The reader is cautioned in the use of these forms to assure that the constraints of continuity 
and the adherence of the index /composition laws are met as required in the derivations. 
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(6) The selection of appropriate initialization functions when structuring fractional 
differential equations will be analogous to, but somewhat more demanding than, the selection of 
constants when structuring ordinary differential equations. Experience working with particular 
equation types and/or types of physical systems will be required of the analyst in the selection of 
initialization functions. 


7. Applications 

The domain of applications and potential applications for the fractional calculus appears to be 
very broad indeed. However, the depth of research and application of the fractional calculus 
appears to be quite limited in most fields, with the exception of rheology (and viscoelasticity in 
particular). Broadly the applications and potential applications are found in diffusion processes, 
electrical science, electrochemistry, rheology, material creep, viscoelasticity, dynamics, and 
controls. Preliminary investigations have also been made into the applicability of the fractional 
calculus to acoustics, boundary layer theory, and turbulence. Further research has been directed 
to the relationships between fractional calculus and fractals. A special (non-exhaustive) set of 
application references has been provided following the general references that direct the reader to 
these areas. It should be noted that not all of the references in this section contain explicit 
reference to the fractional calculus, but are included because either theoretical or experimental 
responses are observed that infer fractional calculus based behavior. 

This section studies a variety of applications and analogues of the operators and combined 
operations of the generalized (fractional) calculus. The intentions are: (1) to demonstrate the 
broad array of uses of the generalized (fractional) calculus, (2) to clearly delineate the effects of 
the initialization function, (3) to contrast generalized versus integer order (calculus) integration 
and differentiation, (4) to demonstrate the generalized zero property, (5) to demonstrate the uses 
of some of the Laplace transform tools (forms), (6) to demonstrate some of the unusual aspects of 
the mathematics, and (7) to demonstrate that the generalization of the zero property, (i.e., eq. 
(5.3.12) may allow the analysis of otherwise insoluble (with the fractional calculus) problems). 
The problems, hopefully, will appeal to both the mathematician and the engineer/scientist. 

7.1 Preliminaries 

Two semi-infinite electrical lines and the op amp (operational amplifier) will be used as key 
building elements for some of the examples that follow. The time domain behavior of these 
elements is summarized here. 
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7.1.1 Half Order Element 


The half order element (semi-infinite lossy line (Hartley and Lorenzo (1998))) is based on the 
one-dimensional diffusion equation 


d v 

~Tt 


= a 


d 2 v 
dx 2 ’ 


(7.1. 1.1) 


which is depicted by a ladder of discrete resistors and capacitors as shown in figure 7-1-1. In 
actuality the resistance and capacitance in this and the inductance in figure 7-1-2 are 
continuously distributed over distance. The terminal characteristics (fig. 7-1-1) of this equation 
(the driving point (impedance) solution) is mathematically described by 


v(/)= rja- 


.- 1/2 


i(t) 


dt 


- 1/2 


+ <Pi (0 


(7.1.1.2) 


or 


v(t) = ryfa c D;' n i(t), where y/(i ,-'/ 2 ,a,c,t) = — 7 =<p,(f) , 

ry/a 


or 

/X 1 d U2 v(t) , , 

i(t) = ~ +<P 2 U) , (7.1.13) 

or 

i(t) = —j= c D', n v(t), where y/(v,y 2 ,a,c,t)= r4a,(p 2 {t). 
r^/a 


■ ■ 
■ 0 • 



V(t) 

Figure 7-1-1. — Half order element. 


Here v(t) and i(r)are the voltage and current, respectively, at the terminal of the element, r is 
the resistance per length of the line, and a is the product of r and c (the capacitance per unit 
length of the line). The y / ' s are the initializations and in both cases are determined by the initial 
state of charge and voltage or current that exists on the infinite array of elements. 
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7.1.2 Zero Order Element 


The second element that will be used is the zero order element (semi-infinite lossless line), 
figure 7-1-2. It has terminal characteristics which are wave equation based and is mathematically 
described (Hartley and Lorenzo (1998) eq. 42) by either of the following equations: 


v(») +«.,(<). (7I-2.D 

or 


i(t)=Jj ,£>(/), where <//(v,0,a,c, 


or 


'(0= i(t)+(p 4 0) . 


or 



(7. 1.2.2 ) 


v(t) = 


c D°i(t), where 


!//(/', 0, a, c,t) = 




v(t) 

Figure 7-1-2. — Zero order element. 


Here £ and c are the inductance and capacitance per unit length of the semi-infinite line. Again, 
it is the initial conditions on the distributed inductors 's) and capacitors (c's) along the infinite 
lines that give rise to the initialization functions (of time). Details of the physics relative to the 
expressions for the <p functions are given in Hartley and Lorenzo (1998). 

7.1.3 Operational Amplifier 


The operational amplifier is a common element in electrical circuits. To analyze it we 
introduce the concepts of impedance and transfer function. For convenience, the circuit of figure 
7-1-3 will be analyzed in a general way here, for use in the later examples. The terms Z, and Z f , 
the input and output impedances (voltage to current ratios), respectively, will represent general 
linear circuit elements together with their initialization functions (therefore are not strictly 
speaking impedances). 
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Then, as Laplace transformed variables 


v,(j)- v g (s) = i,(s)Z i (s) 

input element , 

(7.13.1) 

v tt (s)-v i> (s) = i f (s)Z f (s) 

feedback element , 

(7.13.2) 

o 

11 

1 

1 

grid node . 

(7.1.33) 



Figure 7-1-3. — Operational amplifier circuit. 


These amplifiers are designed so that ideally i = 0 thus v s = 0. Then, after algebraic reduction, 
the resulting (voltage-to-voltage ) transfer function 


y_M = 

v ; (jf) z,(s) 


(7.13.4) 


summarizes the behavior of the device. In the examples that follow we will sometimes use this 
result rather than rewrite the detailed equation set for the elements and node each time. In the 
next five subsections these basic elements will be used to assemble voltage based (computing) 
circuits of both fractional and integer order devices. Several approaches will be used in these 
subsections. 


7.2 Integrators 

7.2.1 Classical ( Integer Order Calculus) Integrator 

An idealized implementation of the classical integrator used in analog computation is shown 
in figure 7-2- 1 . This integrator characterizes integration in the integer order calculus context. The 
defining equations are 

v f (/)-v,(/) = i,R, (7.2.1.1) 
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V„W- v.,W =^, o,-'i / (i) = ^} i,(rVT+ijfi ( (T)rfi 

= 1 1 f/ (T)rfr + ^ = £ j ■> (T)rfr + (v,(c) - v„(c)), (7.2.I.2) 

if = ij and v^,(/) = 0, (7.2.1.3) 

where q(c) is the initial charge on the capacitor. This equation set yields 

v„(r) = -^j v i (r)dr + v <l {c) = —^ e D;\(t), 

c 

with \]/(y i ,-l,a,c,t) = —RCv 0 {c\ (7.2.1.4) 

Then, with RC = 1, this integrator would be represented in the generalized calculus as 
v (f) = _ /), ' v,(/) . The full effect of all past history is then contained in a single number, the 
initialization function, and is the constant - v„(c), compatible with the integer order calculus. 
This is contrasted with the following case. 

c 



Figure 7-2-1. — Classical (order 1) integrator. 


7.2.2 Generalized ( Calculus ) Integrator 


Consider the op amp circuit shown in figure 7-2-2. The input element to the op amp is the 
lossless transmission line (zero order element). Its behavior is described by the wave equation 

d ~ v(x, /) 1 d v(.x, /) ( 720 J 

dt 2 Ic d x 2 

with appropriate boundary conditions. Here the solution at the terminal, equation (7. 1.2.1), will 
be used directly. The defining equations for this integrator are then 

022.1 


NASA/TP— 1 998-2084 1 5 


67 




for the zero order input element, and 


1 1 f 

v , (0 - v o (0 = D;'i, (/) = - J i f dt - v„(c), (7.22.3) 

C 

v s =0. (7.22.4) 



Figure 7-2-2. — Generalized (order 1) integrator. 


Therefore, solving for v„(/) gives 


(*) = - — J Jjv,(/) + p 3 (t) dt+vjc). 


C 


j t 

J J (p 3 (t)dT+ V„(c), 





i (7.2.22 5) 
(7.22.6) 
(7.2.27) 


(7.22.8) 


Taking unity values for the circuit parameters C and c/ £ , the generalized calculus expression 


for the circuit is then 

v„(f)=-.D; l v 1 (r). 


( 7 . 22 . 9 ) 
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This is, of course, the same expression as seen in equation (7.2. 1 .4) for the classical calculus 
based integrator. The important difference is in the values of the respective initialization 
functions. It is further important to note that the detailed evaluation done above to obtain y/ 
would seldom be required in application. For this distributed integrator, the effect of the past 
history is contained not only in the constant v 0 (c), which is the charge on the capacitor, but is 
also carried in the remainder of the y/(t) function, which accounts for the distributed charge 
along the semi-infinite line. It is also observed here that the zero order input element, since it is a 
wave equation, will simply propagate any perturbation in v^/) along the semi-infinite line never 
to be seen again, the only effect being a proportional variation in the current /,(/) . This behavior 
is true for terminal charging, however for the more general case of a side charged line a (an 
additional) time function may be returned to the circuit output, which is dependent on the initial 
voltage distribution on the line. 


7.3 Differentiators 

7.3.1 Classical ( Integer Order Calculus) Differentiator 
An idealized circuit for differentiation is shown in figure 7-3-1. The defining equations are 


= ij i,(T)dT+^ = (73.1.1) 

c < 

v..ir) - v„(l) = i f (t)R . (73.1.2) 

i,(l) = /,(»), v,(l) = 0. (7.3.13) 



This yields 

v„(, )=-«(, (0 = -«c{4(v,(<)-v,.,W)}, (73.1.4) 


Figure 7-3-1. — Classical (order 1) 
differentiator. 


v o {t) = -RC c D]v,(t), t > c . 


(7.3.1.5) 


The initialization function i/r(v f ,l, a. 0, /) is the term — v,_ s (c) , and is normally taken as zero. 


however, it can give rise to a pulse response at time t — c . This classical differentiator 
( RC = l) is contrasted with the case below. 


NASA/TP— 1 998-2084 1 5 


69 




7.3.2 Generalized (Calculus) Differentiator 


The following equation set describes the generalized differentiator as shown in figure 7-3-2. 
1 1 1 

v,(t)-v g (t) = - l .D;'i i (t) = -ji i (T)dr+v i _ i; ( c ).a.3.2.\) : 0 ‘ 


For the zero element 


C 


v g 0 ) - k , 0 ) - cA 0 '/ (0 - ( 0 + 

17 i \ 

J-mi 0, a, c,t ) , (7.3.22) 

V c w y 

f i (r) = * r (r), v„=0 (7.3. 2.3) 


These equations then yield 



V "(') = S | C i'( V ' W ” v »-,(c)) + V'O’/ ,0,fl,c,/)| (7.3.2.4) 


Figure 7-3-2. — Classical (order 1) 
differentiator. 


or 




o,'v,W 


where the initialization function associated with the generalized derivative is given by 


y/ (v n \,a,c,t) = . 


(73.2.5) 


(73.2.6) 


(73.2.7) 


It is noted again that for the terminally charged integer order differentiation, y/ djff = 0 . This is 
proven in the derivation associated with definition of the generalized derivative. As in the 
generalized integrator case, this may also be observed from the physical behavior of the zero 
order element; this is wave equation based and it is well known that disturbances at the terminal 
end will be propagated away from the terminal and are not reflected back because of the infinite 
length of the line. Again, for the more general case of a side charged line a (an additional) time 
function will be returned to the circuit output. 
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7.4 Zero Order Operators (ZOOs) 


7.4.1 Classical (Calculus) ZOO 


An electrical circuit for a voltage-to-voltage ZOO is presented in figure 7-4-1. The defining 
equations are 

v,(,)-v„(r) = R,/, (r). (74X1) Rf 

v„«-v„W = R, /,(/), (7.4.12) 

/,(») = v,(/) = 0. (7.4,13) 

Then from the above equations or from equation 
7. 1 .3.4, with Z, = R f and Z, = /?, , 

v »=-^v,(r)=-^ (7.4.14) 

where \j/(v i ,0,a,c,t) = 0. Clearly this device has no 
memory. Contrast this with the following cases. 



Figure 7-4-1. — Classical (calculus) 
zero order operator. 


7.4.2 Generalized (Calculus) ZOOs 

There are several (many) possible electrical realizations for the zero order operator. Three of 
these are shown in figure 7-4-2, only the last circuit of figure 7-4-2 will be analyzed. Not 
surprisingly, the voltage-to-voltage zero order elements are based on the voltage-to-current zero 
order element described earlier (eq. (7. 1.2.1) and (7. 1.2.2)). The defining equations are 


where 


where 




yr(i f ,0 ,a,c,t) = 


V 


(0-v g (0). 

(7.4.2.1) 


(74.2.2) 

V c i 


(74.2.3) 

r-p 4 (0. 

(7.4.24) 

' / 

v.W=o. 

(74.2.5) 
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Then 


v„(f)=- 



cD?v,(t) 



, oy.it) 



- v s ,0 ,a,c,t 






\t,c 


f ' .JO 


'Cft' 


‘tyXO-.H-VsiO-vAO’ 


and 




£ c 


'/ c * 




The associated initialization is given by 


iff(v i ,0,a,c,t)= (v, - v^,0,^,c,/)+ — y/(i f ,0,a,c,t) 

V C i 


Lc 


= J—' 9> 3 (0 + , t-^» 4 (0. 


C. 


^ c 

c y c / 


(7.4.2.6) 

(7.4.2.7) 

(7.4.2.8) 

(7.4.2.9) 


(7.4.2.10) 


and is the initialization function for the full circuit. Here it is seen that the zero order operator in 
the general case returns the input function v,(f) but also provides the extra time function 
w{v, ,0 ,a,c,t). 



Figure 7-4-2. — Three realizations of the generalized zero order operator. 
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7.5 Half Order Circuits 


7.5.1 Semi-Integrator 


The circuit shown in figure 7-5-1 performs the function of semi-integration of the input 
voltage v, (r). The component equations are 


v,(t)-v g (t) = i i (t)R 

v, (/) - v 0 (t) = r4a c d;'% ( t ) , 


(7.5.1. 1) 

(7.5.1.2) 


where y/(i-\/2,a,cj) = —^=(p ] (t), (7.5.1.3) 

rand a are as defined for equation (7.1. 1.3), 
and v^=0. (7.5.1.4) 

Solving for v (( (/) 

v„(») = -rVS ,.D,- ,,! 4v,(r)j, then (7515) 




where 



Figure 7-5-1. — Semi-integrator. 

(7.5.1.6) 


i i/(v j ,-/ 2 ,a,c,t)= Ry/(i f -/ 2 ,a,c,t) = —j=(p ,(?). (7.5. 1.7) 


Here equation (7.5. 1.6), with the leading coefficient specialized to one, is the basis of the semi- 
integrator computing element. Note, the equivalent (uninitialized) impedance form may also be 
readily found from equation (7. 1.3.4), because Z f = r*J a / s l/_ and, Z, = /?, this is given as 


v„(s) _ _rja_ 
v,W" Rs m ■ 


(7.5.1.8) 
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7.5.2 Semi-Differentiator 


The circuit of figure 7-5-2 yields an output voltage proportional to the semi-derivative of the 
input voltage. For this circuit, 

'■(0 = ^ ( D,'' ! (v,W-v t (/)), {1523) 

where rand a are as defined for equation (7. 1.1.3), 


v{ v i -v H ,y 2 ,a,c,t) = r4a(p 2 {t), (1.522) 

v g (t)-v„(t) = i f (t)R, and (7.5.23) 


Then, 


i,(t)=i f (tl v K (t) = 0. (75.2.4) 

(7.5.25) 
(75.2.6) 


v.W— «/,(»)—- 

ryja 

mdy/ (v n y 2 , a, c,t) = yr(v ( . - v,,, l / 2 ,a,c,t) =rVa (p 2 (t). 


Here equation (7.5.25), with the leading coefficient specialized to one, is the basis of the semi- 
differentiator computing element. These theoretical circuits then can be used as the basis of an 
infinite variety of combination circuits. A few of these will be examined to illustrate some of the 
issues of the earlier theoretical development. 



Figure 7-5-2. — Semi-differentiator. 
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7.6 Combination Circuits 


7.6.1 Tandem Semi-Integrators 


Consider the circuit shown in figure 7-6- 1 . For this circuit, 

r,./«7 

v a (0 =- » 

K \ 

with y/ ] (v, }4, a, c,t), and 

r 7 -lcc^ 

with \]/ 2 (v 2 -/ 2 ,a,c,t) . 


Then 


r,r, Ja,a 2 , , 

v 3 (0=- l ^t; — .•*>>,(/), 




R 

with i/(v p — l,a,c,/) = f d~' n \lf x {y v — X A ,a,c,i) — j=\ ff 2 (v 2 ,~ 'A,a,c,t) 

r 2 \ ®2 



Figure 7-6-1 . — Integrator composed of tandem semi-integrators. 


(7.611) 

(7.6.1.2) 

(7.6.1.3) 

(7.6.1.4) 

(7.6.1.5) 

(7.6.1.6) 
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7.6.2 Inverse Operations 


Consider the circuit shown in figure 7-6-2. For /?, = R 2 = r, = r 2 = = ^/a7 = 1, the 


solution is 


>M =-o! , 2 (-d;'\(')) - , 


as expected. Note however, that v 3 (t) = v,(/) only under the condition that 
y/(v j ,0,a,c,t) = C d'/ 2 y/(v { ,-y 2 ,a,c,t)-y/(v 2 ,/ 2 ,a,c,t) = 0 


where y/(v,,0,a,c\?)is the initialization function for the combined operation (i.e., c £>,° v, (r) ). Under 
terminal charging, of course this is satisfied. However, under side charging conditions it is seen 
that v 3 (r) = v,(r) + ,0,a,c,t) and the zero property is not (necessarily) satisfied as indicated 

previously. It should be clear from this example however that the ability to readily analyze such 
systems is a perfectly acceptable (indeed desirable) state of affairs. 

The discussion will now move away from circuits analogous to generalized calculus 
operations and consider applications from other areas of engineering and science. 




Figure 7-6-2. — Inverse operations using semi-operators. 
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7.7 Dynamic Thermocouple 


It is desired to determine the dynamic (frequency) response of a thermocouple (fig. 7-7-1) 
that is designed to achieve rapid response. The thermocouple consists of two dissimilar metals 
with a common junction point. To achieve a high level of dynamic response, the mass of the 
junction and the diameter of the wire are minimized. Because the wires are long and insulated 
they will be treated as semi-infinite (heat) conductors. The analysis will determine the transfer 
function relating the junction temperature to the free stream gas temperature T h (s ) / 7j,(s) . For 
the semi-infinite conductors the conducted heat rate Q(t) is given by 


< 2,(0 = 



d U2 t, 


h ’ 


(7.7.1) 


where k is the thermal conductivity and a is the 
thermal diffusivity. For the transfer function the 
effects of initialization are not required, therefore, 

allyr(0's are zero. Thus the following equations 
describe the time domain behavior: 


Q t (t) = hA(T t (,)-T h (t)), (7.7.2) 

T„(t) = — ,D;' (Q,(i) - 0 (f) - Q, (/», (7.7,3) 


wc v 


Q,(l) = -j== J>', l2 T h U \ , and 


Jcc. , 


(7.7.4) 

(7.7.5) 


where h A is the product of the convection heat 
transfer coefficient and the surface area and wc v 
is the product of the junction mass and the specific 
heat of the material. Taking the Laplace transform 
of these equations, eliminating the <2's, and 
collecting like terms yields 


wc v s + h A + 




hAT y {s), (7.7.6) 


i < 2,(0 



Figure 7-7-1. — Dynamic thermocouple. 




Figure 7-7-2. — Frequency response of 
dynamic thermocouple. 
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which gives the transfer function as 


r»M i 

wc v ) 1 


s' n +l 

\hA ) hA 



(7.7.7) 


The magnitude and the phase angle are readily determined by letting s = ion and noting that 
V / = cosk/ 4 + i sin 7r/4 is the principle value. Figure 7-7-2 shows the transfer function in the 
form of a Bode plot for 


wc\. 

= 0.005 and 
hA 



= 5.0 . 


The value of the generalized (fractional) calculus is clearly demonstrated in this application; 
conventional approaches here would require the solution of two simultaneous partial differential 
equations with an ordinary differential equation. While this is straightforward, the effort and 
required attention to detail would be much greater. The response shows two distinct asymptotes; 
in the mid- frequency range a slope of - 10 db/decade corresponds to the 1 / ,v l/2 behavior and a 
slope of -20 db/decade for frequencies above 10 6 radians /second corresponds to 1 / ,v behavior. 


7.8 Electrolytic Cell 


The electrolytic cell is well known to exhibit fractional behavior. Typical analysis is based on 
impedance methods which ignore initialization effects. This example analyzes a simple model of an 
electrolytic cell (battery) (fig. 7-8-1). The fractional element is a half order system representing the 
electrode-electrolyte interface. This is a diffusion process and is known in that field as the Warburg 
impedance or constant phase element (see for example, Bard and Faulkner (1980)). Here we have 
included an initialization along with the basic properties of the Warburg element. Two phases of 
behavior are considered: ( 1) a charging/relaxation phase and (2) a load drawing or usage phase. The 
charging/relaxation phase takes place between t = a = 0 and t = c , with actual current flow 
(electrical charging) only occurring for a = 0 < t < b . Later, in the load drawing phase analysis the 
desired charging response may be time shifted to set c = 0 , if that is desired. The end point c, 

(c > b ), in the charging analysis will be kept variable and may be quickly obtained by setting t = c 
in the solution equations of this section. 
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7.8.1 Charging Phase 


The time domain equations for the charging/relaxation phase (fig. 7-8- 1(a)) are 


v i(0" v 2(0=^oA"'*v(0 

(7.8.1. 1) 

= ^Jl(^ / + V .-2( 0 )’ 

^ 0 

v 2 (t) — V 3 (f ) = ’ 

(7.8. 1.2) 

*V (f ) = * J, W + hyM, and 

(7.8. 1.3) 

■VW=i„A ,,! (v,(>)-v,(rt). 

(7.8. 1.4) 

Taking the Laplace transform yields 


( \ it 1 • i u 

v,(s) — V,(s) — l ( \S) + 

Cs s 

(7.8.1.5) 

v 2 (s)-v ? (s) = i R {s)R , 

(7.8. 1.6) 

( (•*) = (^) “I" ( 5 ) * and 

(7.8. 1.7) 


hv (‘0 ~~ g ' (' 2 (“') ^3 (^)) (^2 V 3»/2 ,0, C, i)) ^ g j 


WAR 


V ^H 


V)(0 


v->(0 


iw<0 


ic'(t) 


HV\H 


<R«) 


(a) Battery - Charging phase circuit 


rlH 


icd) 


WAR 


v 2 (t > i w (t) 




v^(t) 


iR(t) 






(b) Battery - Load usage phase circuit 




-1 t 


a=0 b c 

(c) Charging current profile 


... , . . ... . / x . / -v , 1 \ Figure 7-8-1. — Electrolytic cell circuit. 

Now eliminating the internal variables *», i R {s) and v 2 (s) (a) Charging phase . (b) Load usage phase . 

and rearranging, yields the following forms Charging curTent P rof,le - 


((•v) = 


''(i ■' 


-v ,/2 + i) 


v + _L + 

A T CB A ^ RC 


r( v i( 5 ) - v j( 5 )) 


+ 


^sy/(v 2 - v 3> K,0, £,■*)- V|- 2 ( 0 )(l 5 ' /2 + «) 


, 1 J/2 . 1 

5 + ^5 +- 


/?c 


or v,(.v) - v 3 (j) = 


5 + ^« 1/2 +^f 




■ s ' l/ ' + r) 


{(s)- 


j I f/(s) V,_ 2 (0) 


1 _l/2 . i 
B A ^ ff 


(7.8.1.9) 

(7.8.1.10) 


NASA/TP— 1998-208415 


79 




These forms can be used either for transfer function analysis or to determine time responses 
including the effect of initialization. The following transform pairs apply to the forced term of 
equation (7.8.1.10) 



(7.8.1.1 1) 


(7.8.1.12) 


(7.8.1.13) 


During the initialization period, 0 = a < t < c , both of the initializations are taken as zero. Then, 
for example, with impulsive charging, i c (s) = 1 , and 


v', 






(7.8.1.14) 


For general /(/), the convolution theorem is applied yielding 



(7.8.1.15) 


This is the solution (response of overall cell voltage to charging current) for the charging phase 
behavior for a general / (t) . If the current is taken as a constant during the electrical charging period, 
/ (/) = / , for 0 = a < t < b, then, / (/ - r) = /,(«(? - t)-u{t -T-b))(see fig. 7-8- 1(c)) and this 
equation becomes 


v,W-v,W = 


J 






t> 0. (7.8.1.16) 
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The voltage response v,0) - v,(/) of the battery during the electrical charging period 0 < t < b = 1 and 
during the relaxation period t > 1 are shown in figure 7-8-2. Also shown is the voltage response of the 
Warburg element v,(/) - v,(?) (determined in subsection 7.8.2). The battery parameters are taken to be 

R = B = C= 1.0. 



Time 


Figure 7-8-2. — Voltage response with time for electrolytic 
cell-charging and relaxation phases. 
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7.8.2 Initialization Function 


The above response (eq. (7.8.1.16) or (7.8.1.18)) then can be the basis for a “good” 
mathematical form for an initialization function for systems (cells) of this type with constant current 
charging. It of course includes both the effects of the Warburg element and the capacitor. These 
may be separated as follows. Since v,_ 2 (0) must be zero, the voltage across the Warburg element 
v 2 (r) — v 3 (/), is determined as 

v 2 (0 - v 3 (f) = (v, (t) - v 3 (/)) - (v, (t) - v 2 (tj) = (v, (/) - v 3 (/))-^J i c dx , t > 0 , 

= (v, ( t ) - v 3 (f))-^j /,(t< t) - 1{ r - b))dr , 

which using equation (7.8. 1.16), yields 


=1® 


\ TIT R l 


b4t 

R 


-Idr~ 


-b 

1 B (Jfr 1 

(Bji) 

H 

0 

r - e ef fc 

4kt a 

M 


I dr, t> 0. 


(7.8.21) 


This is the voltage across the Warburg element during the charging period. Note, this 
equation also applies for t > c with an open external circuit. The input function to (argument of) 
the semi-derivative of the Warburg impedance is (v 2 (f)- v 3 (/)) . If it is assumed that the Warburg 
impedance is terminally charged, then the initialization function for the Warburg element may be 
determined from equation (6.3.3), namely, 
d m 

il/(f,u,a,c,t) = -j^-l(f(f,-p,a,c,t) + y/(h,m,a,c,t). (6.3.3) 

Because terminal charging is assumed = 0, thus, 

\ d ^ 

^^«(v' 2 -v 3 ,)<,a,c,/) = -p^— J(?- t) K ''(v 2 (t)-v 3 (t))^t, t>c. (7.8.22) 


This then is a “good” initialization of the Warburg element in the context of this electrical 
configuration and charging profile. This initialization applies for t > c> b>a = 0. In its 
application, if it is desired to initialize with c = 0 , then it would be necessary to time shift this 
result appropriately, that is, replace t by t - c in equation (7.8.2. 1). Shown here is the formal 
approach to determination of the initialization function. In the next section, material creep, an 
alternate approach, which is sometimes useful, will be shown. 
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7.8.3 Load Drawing Phase 


Consider now the load phase behavior of the electrolytic cell as shown in figure 7-8- 1(b). The 
defining equations are amended to reflect the added external load resistance. 


v, (f ) — ,D;'i (t) = -^;]i,(t)dt + v,_,(c), (7.83.1) 

v, (<)-»,(<) =;„(')«, OS32) 

K (0 = h (0 + hr (0 = hi. (0» (7.83.3) 

v 7 ,(t)-v,(t) = i RL {t) R L ’ and (7.83.4) 

UO'-J^rMO-v^)). (7.83.5) 

Taking the Laplace transform yields 

v,W-v 2 W = ^i,W+— . (7.83.6) 

Cs s 

v 2 {s)-v 3 (s) = i R (s)R , (7.83.7) 

i t (s) = i R (s)+iw(s) = /«-(*), (7.83.8) 

v 3 (j)-v l (5) = / 1M .(j)/J t , and (7.83.9) 

^ (s 1 ' 2 ( v 2 (s) - v, (.?)) +y WAR {v 2 - V, , >< ,0, c, ,v)) . (7.8.3.10) 


Then from the charging analysis, equation (7.8. 1 .10), with i c (s) — i RL (s ) — (v 3 (s) v, ( 5 ))/ R L , and 
taking a = 0, and c = 1 , 


/ 

R l s 

<B + Rj 

\ 

V,_ 2 (l) VwAr{v2-Vi,'/2,0Xs) 

0 ^ 

R , i- 

B 

( R,\ 

1 + -^ 

R 

1 X 1 

5 4 S z - + \~~ 

CB RC J 

S 

B 

V 

V l' 

, B + R, 

/ 


This response is seen to have no forced term and two initialization driven terms. These are 
the initialization effect of the capacitor, i.e., the v,_ 2 (l) term, and the historic effect of the 
Warburg impedance as reflected in the y/ WAR term. It is readily analyzed using the following 
transformation, let p-s ^ => p 2 = s, = s ^ , then 


v,(p)-v 3 (/>) = 


/ 

P~ 

\ 

B ) 

p+ — 

R) 

4 

/ j \ 

Vi ,(l) g^(v 2 -V 3 ,/ 2 ,0,l ,/> 2 ) 

„3 * 

f R I 

1 + — 
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L*. 1 d+— 

P 2 ' z + ! 

B R ) 

l R l 

) p r l c p rr l c 


. (7.8.3.12) 


NASA/TP— 1 998-2084 1 5 


83 



From this point, an appropriate substitution or approximation is made for y/(p 2 ) , and 
conventional analysis tools (factoring, partial fractions, etc.) are applied to evaluate the p 
domain response. This response may be transformed back to the s domain and subsequently 
inverse Laplace transformed back to the time domain if desired. This is expedited by the use of 
special transforms, which will be presented in a later paper. 

7.9 Material Creep 

7.9.1 Creep Strain 


Creep behavior of materials may be described and analyzed using the methods of fractional 
calculus. A stress-strain law for certain viscoelastic materials has been suggested by G.W. Scott- 
Blair (see Bland (1960) p. 5). It neglects both the instantaneous elastic and the long term viscous 
flow effects. The relationship in terms of this paper is given by 


e(t) = j r D r^(t) 

or, when the specimen is new, no initialization is required hence, 

e(t)=^.drm 


(7.91.1) 


(7.9.1.2) 


Here e and a are the strain and stress respectively. K and v are constants associated with 
the material, v = 0 for an elastic solid, and v = 1 for a viscous liquid. 

Consider a uniform stress load to be applied between times t = 0 and t = d , on a new 
specimen (no initialization) then. 


(7.9.1.3) 


where H is the unit step function, and k is the magnitude of the load. 
Then, 


e(t) = — j (t - t [tf ( r ) - H( r - d)]d r , 

A1 (V)' 


f / 

J(/ - t) v_1 dr H(t ) - J(r - r) v_1 drH{t- d) 


*T(v)[_ 

f—r[t v H{t)-(t-d) v H(t-dj\, 


(7.9.1.4) 
and (7.9.15) 

v , + 1) I ‘ '* v ‘ / v ‘ " v ‘ **/i» (7.9.1.6) 

defines the strain (creep) response for the prescribed loading. The effect of v is shown in figure 

. KT(v + l) 


7-9- 1 , where normalized strain (creep) is 


K 


£ (t) . The response £ ( t ) for t < d = 1, is 


seen to be essentially the creep function. The response for t > d = 1 is the stress relaxation period. 
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Figure 7-9-1 . — Creep (strain) versus time for various v. [v = v in figure] 


7.9.2 Creep Initialization 

For creep applications of this type the initialization function associated with constant past 
loading may be readily inferred from the solution for equation (7.9. 1 .6), as follows. 

Consider any problem initialized at point c, then, 

,D, -'/(<)=., 0;’/(>). '><■. (7.9.2.1) 

XVM+v(/.-«-«.*.0=.^/M=xv(0. 092 . 2 > 

thus y/ may be expressed as 

y/(f -q,a,c,t)= a d~ q f(t)- c d; q f{t) , t>c . (7.9.23) 

It is now desired to obtain an initialization function for the creep problem initialized at 
t = c = d where d is as used in equation (7.9. 1.6). Applying this result, equation (7.9.23), to the 
creep problem gives 

c D~ y G{t)=d~ v c{t)+y/(o-v,a,c,t), t>c. (79.2.4) 

Then taking c — d and since o{t) - 0 for t > d, 

,. d;'o{t)= d d;'(j{t) = 0 , t>c-d, (79.2.5) 

and K times the response equation (79.1.6) is the initialization t/r(r)thus, 

= t>c = d. C 2.92.6 ) 
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This initialization function may now (with proper time shifting) be used for specimens which 
have experienced a creep history. 


7.9.3 Relaxation Function 


The stress (creep) relaxation function ^(t)is defined as the stress required to produce a strain 
H(t), excluding terms that are initially infinite or do not tend to zero as t becomes infinite 
(Bland ( 1960) p. 38). This is easily determined in the fractional setting. For a new material, the 

stress is given by 

* 

X(t) = o(t) =K,d;e (r) , 

(7.9.3.1) 

■**ni-v)! ( '" T) w<lVr ’ 

( 7.93.2) 


(7.93.3) 


yielding as a result the creep relaxation function for this formulation. 


7.10 Fractional Order Tracking Filter 


The availability of fractional differintegrals allows new freedoms in many areas, tracking 
filters is but one example. Here a noisy signal x(t) is to be filtered to yield the filtered signal 
y(t ) . Only a single active element will be considered here, then the filter is described by 

( 0 D?+a)y(t) = ax(t) (7.10.1) 


or 

0 d^y(t) + y/(y,q,0,c,t)+a y(t) = ax(t) . (7.10.2) 

Taking the Laplace transform yields 

s q y{s) + y/{s) + ay(s) = a x(s) (7.10.3) 


or 


y(s) 


a 

s' 1 + a 


44 - 


vU) 

s v + a 


(7.10.4) 


While the y/ function may allow the possibility of “precharging” the filter to minimize the 
filter lag, the initialization function in this application will be discarded. The transfer function of 
the filter then is given by 


44 _ £ 

jc(i) s‘‘+a 


(7.10.5) 
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Then the behavior of the filter can be considered by studying the response to a unit step in x(t) thus, 

a 


yU) = 


+ a) 


(7.10.5) 


The full range 0 < q < 1 is of course available, and requires a series solution for the inverse 
transform. This will not be detailed here, but time responses for a step input for *(/) showing the 
effects of varying both a and q are presented in figure 7-10-1. The inverse transform for the 
special case q = 'A is readily shown as 


y(t) =\ - e‘ r> Erfc(a t' A ) . 


(7.10.6) 


The frequency response for the tracking filter (eq. (7.10.5)) for various q with a = 1 , is 
presented in figure 7-10-2. 




Figure 7-10-1 . — Fractional order tracking filter response to unit step input, 
varying q (a = 1 ) and varying a (q = 0.5). 
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Figure 7- 1 0-2. — Fractional order tracking filter frequency response, varying q and a - 1.0. 


7.11 Viscoelastic Damped Vibration 

The addition of damping is important in many applications for stability augmentation. One 
such application in the aerospace field is damping augmentation for gas turbine fan blades. The 
vibration modes in this case are typically bending, torsion, or combinations. Here it is desired to 
determine the analogous linear vibration mode behavior. The viscoelastic damping may be added 
in different ways. The interest here is in determining the transfer function (.v domain response) in 
general terms, then, by specialization of the constants, various cases may be examined. 
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Since the transfer function requires the initializations to be zero, the equations describing the 
generalized system (fig. 7-1 1-1) are 


i 

* 

1 

0 

^3 

c 

1 

ii 

(7.1 1.1) 

H 

1 

i 

II 

tc 

(7.1 1.2) 

^2 ” ^2 0 df (X 0 X m2 ), 

(7.1 1.3) 

^2 = ^3 od/' (-* m 2 

(7.11.4) 

F ] + F 2 = rn 0 dfx 0 . 

(7.11.5) 

Where the k’ s are damping coefficients and the spring 
constant, and m,F, and x are the mass, force and 


displacements respectively. 



Figure 7-1 1-1. — Spring-mass-visco 
damped dynamic system. 


Taking the Laplace transform of these equations and eliminating x m] and x m2 yields 


A 



( \ 

-ky 1 

(x 0 -x,) and 

F 2 (s) = 

- ky 2 

l + ^-s* 1 
k ) 

l + h. s M 
k ? J 


(7.1 1.6) 


Then after eliminating the forces and some algebra the generalized transfer function is obtained as 


-r n (s) 


_L + 1 — x ‘i t+, i i+ ‘i i ‘ 

[■ k k k 

22 22 2! - — - : .(7.1 1.7) 

m s 2+<jl + m + m s 2+,,2 + m s l+<i>ui2 + JL a .</!+^ + _L v '/ u ‘/" _| s <lW _| 5 «I+«2+i/J 

k { k 2 kk 2 k t k 2 kk 2 k 2 k ] k 


Now various special cases may be examined by appropriate choices for k,k n and the q s . For 
example, to allow damper 3 to represent a conventional (dashpot) damper, let q3 = 1 and select 
the appropriate k y Here, for simplicity, the glandg2 dampers will be eliminated (made into 
springs of infinite stiffness), that is, q\ — q2 = 0, and k t = k 2 = », and the q3 damper will be 
considered to be viscoelastic (i.e., fractional value for q3)\ then the transfer function for figure 7- 
11-1 is given as 

k k 

*3 r/3 . * 

( \ S + 

foW m m (7118) 

* 1 
m m 
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Phase (degrees) 20*log10(|Xo(s)/Xi(s)| Phase (degrees) 


The frequency response may now be evaluated by letting 5 = ia> and determining the magnitude 
and phase angle (polar form argument) using equation (7. 1 1 .8). The results for kjm = 1 and 
various values of k^/m and q3 are shown in figures 7-1 1-2 and 7-1 1-3. 



Frequency (rads/s) 


Figure 7-1 1-2. — Mass-spring-viscoelastic damper transfer function, 
varying q3 the viscoelastic coefficient from 0 to 2 by steps of 0.2. 



Frequency (rads/s) 



Figure 7-1 1-3. — Mass-spring-viscoelastic damper transfer function 
K3/m varies from 0.25 to 3 by steps of 0.25, q3 = 0.5. 
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8. Variable Structure (Order) Differintegral 


8.1 Variable Order Integration 

8.1.1 Research Issues 


This section presents some research questions which have been enabled by the fractional 
calculus. Consider the fractional differential equation 

c D?y(t) = f(t) (811.1) 

and the inferred integral equation 


c Dff(t) = y(t). 


(811.2) 


Since q in the fractional calculus can take on any real (or complex) value, the question is asked, 
What is a desirable definition for the fractional integral when q is allowed to vary either with 
t or yl More specifically, What is an appropriate definition for 

f D f ' ,( '- v) /(0 ? (8113) 


8. 1.2 Motivation 

There is considerable potential physical motivation toward the creation and implementation 
of such a concept. A few possibilities will be mentioned here. 

From the field of viscoelasticity (Bland (I960)), the effect of temperature on the small 
amplitude creep behavior (force/extension) of certain materials is to change the characteristics 
from elastic (spring-like, q = 0 ) to viscoelastic or viscous (damper-like, q = -1 ). This relates to 
the expression 

c Df F(t)=k[x a {t)-x h {t)), (812.1) 

where F is force and x is displacement. Experience (and experiments), typically, is based on fixed 
temperature, but real applications may well require a time varying temperature to be analyzed. 

Glockle and Nonnenmacher (1995) studied the relaxation processes and reaction kinetics of 
proteins which are described by fractional differential equations of order (} . The order was found 
to have a temperature dependence. (It is interesting to note that this work also alludes to 
difficulties which appear to be related to initialization). 

Smit, W. and deVries (1970) studied the stress-strain behavior of viscoelastic materials 
(textile fibers) with fractional order differential equations of order a , with 1 > a > 0. They 
show, based on related experiments, a to be dependent on strain level. 

Polymer linear viscoelastic stress relaxation was studied by Bagley (1991). This process is 
described by fractional differential equations of order/? for a given fixed temperature. The paper 
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shows a clear dependence of /? on temperature for polyisobutlene and correlates fractional model 
and experiment. Further, it is indicated that p order fractal time processes lead to p order 
fractional derivative consitutive laws. 

From the field of damage modeling, as the damage accumulates (with time) in a structure the 
nonlinear stress / strain behavior changes. It may be that this is better described with variable 
order calculus. 

Finally, the behavior of some diffusion processes in response to temperature changes may be 
better described using variable order elements than time varying coefficients. 

8.1.3 Analysis 


For present purposes, attention will be focused only on variations of q with t, that is, 
q = q(t) . Further, initialization will not be considered to start, and the Riemann-Liouville 
definition will be taken as the basis (see Appendix C for Griinwald variable structure basis 
discussion). For simplicity in these considerations c = Oand y/(f,-q(t),a,c,t ) = 0. The most 
general form then under consideration is 

( 

0 

where the notation q e (t, r) and q R (t, r) have been used to indicate that the arguments of the q ' s of the 
exponent and the gamma function may be different. The simplest arguments to be considered are 
q(t, r) — » q(t), q(t, r) — > q(r) and q(t, r) — > q(t - t) . For the two occurrences of q in the defining 
integral (eq. (8. 1.3.1)), this yields 9 permutations to be studied. To quickly screen these possibilities and 
more, computer simulations of the defining integrals (with a variety of test functions, /(?)' s ) were 
used to quickly eliminate potential definitions with unsatisfactory properties. The desirable properties, 
of course, are those contained in the criteria examined earlier. Clearly backward compatibility, and the 
zero properties are not issues for definitions which are special cases of equation (8. 1.3.1). This leaves 
linearity and the index law (composition) to be considered. 










-f(r)dT, 


( 8 . 13 . 1 ) 


8.1.4 Criteria 


n , 


Linearity will be considered first using the general form, equation (8.1.3. 1). Then, 
0 D; ll{,) (a f(t)+ bg(t ) ) = ] — -r - {a f(z)+ b g{z)}dt 

' f ( t _ T V'd' 0-> ' ( t _ 

= a 1 11 7 \\ f( r ) dz+b i r 7 vT g ^ dr ' 

0 r(^(t,r)) J 0 r(q A ,(r,r)) 

=a 0 D;^f(t)+b 0 D;^g{t). 

Therefore linearity is satisfied for all arguments of q in the defining integral. 


(8.1.4.1) 

(8.1.4.2) 

(8.1.4.3) 
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The computer simulation showed that the interesting behavior occurred when 
q e (t,x) = q K (t,t). Of these three cases only the case q(t, t)-» q{t - t) provided adherence to the 
index law (for selected functions); but it did not satisfy composition, namely, under the definition 

(8144) 

it was inferred that 


f(t)= 0 d; vU \d;< {,) f{t)± 0 d;« {,) - vU} f(t). 


(8.1.4.5) 


It is observed that this definition is a convolutional form. The convolution character of the 
generalized integral will now be used in proving the above property; indeed the property is fully 
equivalent to the associative law for convolution integrals (see Mikuzinski (1967)). 

The index law will now be proved. The LHS of equation (8. 1.4.5), under the definition of 
equation (8. 1.4.4), is 


oD, 


vU) 


o D, 


->•(/ 


/ 


/(<) = J 


0 


r(«(< - r)) 


J „ r(v(r-j)) 


f(s)dsd T. 


( 8 . 1 . 4 . 6 ) 


Replacing the outer integral with its convolution equivalent yields 


</(r)-l 


r (r-T-v) 


v(r-T-s)-\ 


|r(«W)! r(v(<-r-s)) 


■f{s)dsdx ; (8.1.4.7) 


letting s — 0) — z => ds - d(0, 
r T q{zy ' r (t - 

= (8. 1.4.8) 


J „r( 9 (T))J r(v(r - o)) 


=Ji 


T g(rH {t - (o)' 

ii r{ q (r))r{v(t-(o)) 




f{(0-x)d(Odr. (8.1.4.9) 


The double integral is extended to the triangular region 
as shown in figure 8-1-1, yielding 


>} r( ? (r))r(v(i-®)) 


f((o-x)d(o dr . 


(8.1.4.10) 


T 




Figure 8-1-1 . — Reference 
triangular region. 


Now the double integral is converted back to an iterated integral with the order of integration 
reversed, giving 


! r(v(r-®)) 


! r(«(0) 


f((0-z)dT d(0. 


(8.1.4.1 1) 


NAS A7TP— 1 998-2084 1 5 


93 



(8.1.4.12) 


The inner integral is now replaced with its convolutional equivalent 


K/-Q?) v(f - MH 

| r(v(/-©)) 


r (ft) - -r) ?(<a ~ f) ~ 1 

I r(^(o> - r» 


f(t)dz d(0. 


= 0 D, 


-*'(/) r\-<K0 




(8.1.4.13) 


which completes the proof. At this time the definition of equation (8. 1.4.4) appears to be the 
most satisfactory definition in terms of these properties and its strong relationship to convolution 
theory. It is noted explicitly that composition is not satisfied, that is, that 


, o D~ v(,) f (f ) = 0 


(8.1.4.14) 


However, the lack of this computational convenience does not diminish the potential usefulness 
of this concept for application in the physical sciences. It would however, considerably increase 
the analytical difficulty. 

It should be noted that the definition 

. j ^^ r - f(-c)dz q(t)> 0 (8.1.4.15) 

can be shown to have a weak form of composition, namely, 

0 d,-" (,) 0 D; k f(t)= 0 D; {,l( ' >+k] f(t)* () D; k 0 D; u{,] f(t), u{t\k > o. (8.1.4.16) 
Also the definition 

• /(». «(«)> 0. (8X4.17) 

can be shown to have the following weak form of composition, as 

„D;* 0 D;* l /(»KA' ( ** v<, V«*oiT ( ' > . k,v(,)>0. (8.1.4.18) 


The weak composition results indicated here have been validated by computer numerical 
simulation for a variety of functions. 

The question of which is the most desirable definition for variable order integration, of 
course, is still an open issue. Based on the preliminary investigation of this study, the most 
compelling definition is given by 



r(<7(r-r)) 


f(r)dt , 


q{t)> 0, (8.1.4.19) 


because of adherence to the index law and because of the convolutional form which makes 
available all of the results of the associated theory. It may be possible that a slightly more 
complicated form of composition can be shown than that considered in the inequality (8. 1.4.5). 
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Several candidate forms have been considered without success. It is noted that different physical 
processes may effectively use different definitions and that all three of the forms, equations 
(8. 1.4.4), (8.1.4.15), and (8.1.4.17), may prove useful. Finally, it is further noted that a set of 
differintegrals parallel to those based on the Riemann-Liouville fractional integral (equation 
(8. 1.3.1)) may be formed based on the Griinwald definition. A brief consideration of the variable 
structure differintegral based on the Griinwald definition is given in Appendix C. 


8.1.5 Laplace Transform 


The derivation of the Laplace transform of the variable structure integral follows that for the 
fixed order case exactly since the convolution theorem can be applied. Then, considering the 
uninitialized case 

'r(f- t)* mK 


L{oD,-« (,) /(r)} = J g~”l J /( T ) jT dt . q(t)>0,t>0, (81-5.1) 


then as before 


' i 

L{h(t)* g(r)} = H(s)G(s) = l\ h(r)g{t - z)dr 


(8.1.5.2) 




Now taking h(t) = f(t) and g(t) = ^ ^ the convolution theorem yields 

f(s)gU) = ■ 


(8.1.5.3) 


The variable structure differintegral allows the introduction of a new transfer function 
concept. The conventional transfer (see fig. 8-1-2) function relates the Laplace transform of the 
output to the transform of the input by 


TF l = 




q(l) 


i{/(<)} L- 

[/(' 

)} 

Y 

II 

l r («W)J 



(8.1.5.4) 


/(<) 





Figure 8-1-2. — Block diagram for transfer functions 
of variable structure integral. 
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Since q is now a variable it also may be thought of as an input, and a new transfer function may 
defined as 

(8. 1.5.5) 

For the considered definitions, the process from q(t) to y(r) has not been shown to be linear. 
That is 

o D -v.<.H,2(r> / (,) = o o D~ q2(,) f{t) 

has not been shown. Thus, the meaning and utility of equation (8. 1.5. 5), for these definitions 
requires further consideration. 

The relationship of the two transfer functions may be determined as follows. Consider 
/(/) and q(t) to be related by g(t), where 

t 

<?(') = | - r )dr then L{q{t)} = L{/(/)}L{g(r)} (8.1.5.6) 

0 

by the convolution theorem. Then 

_ L{y(f)} _ L{f(,)}TF, 

1 L{ q (,)} L{/(r)}L{«(r)} 

8.2 Variable Order Differentiation 

If the form of equation (8. 1.4.4) is assumed, then the parallel definition for the variable 
structure derivative might be given (considered) as 

0 Df' ] f(th 0 D? 0 D; p(,) f(t), q(t)> 0, (8.2.1) 

where q(t) = m - p(t), and m is a positive integer. If it is assumed that q(t ) is always positive, 
then m could be taken as the least integer greater than q(t ) . However, since composition does 

not hold it is not at all clear that this is a reasonable definition. Matters are further complicated 
by the fact that it may be desirable to allow q(t) to range over both positive and negative values. 
This places a “seam” at q = 0, which may make any approach based on the Riemann-Liouville 
definition implausible and perhaps require an approach based on the Griinwald definition. 
Consideration of the variable structure differintegral based on the Griinwald definition is 
presented in Appendix C. 
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9. Summary 


The primary goal of the work discussed in this paper has been to make practical the 
application of the fractional calculus to the problems of science and engineering. This has 
required a new (modified) set of definitions for the fractional calculus. This formulation includes 
newly developed initialization functions that bring in the effects of the past for the system being 
analyzed. For fractional integrals and derivatives with respect to time, the initialization functions 
are time functions, which are a generalization of the initializing constants required for ordinary 
differential equations. 

The following was accomplished: 

• Presented several conceptualizations for the fractional differintegrals, both analytically 
based and physically based. Conceptualizations of the fractional differintegrals in terms 
of weighted time delays of the argument f(t) (conveyor analogy) have been developed. 
Geometric interpretation of the generalized differintegral have been developed. 

• Developed a formal initialization for the fractional integral and derivative which is 
incorporated directly into the definitions. 

• Broadened the definition of the integer order derivative to include initialization. 

• Proved that the Ross criteria for a fractional calculus, namely, that backward 
compatibility, zero order property, linearity, and composition hold under terminal 
charging conditions. Provided inter-relationship constraints between initialization 
functions required for satisfaction of the criteria under side charging conditions. 

• Eliminated the requirement that the function and all of its derivatives be zero at the 
starting point, t = c, by the introduction of the initialization function. 

• Developed Laplace transform expressions for the generalized fractional integrals and 
derivatives (with initialization included). 

• Developed Laplace transforms for integer order and fractional order decompositions of 
both the generalized fractional integral and derivative. 

• Determined the effect of the starting point on the Laplace transforms of the fractional 
integral and derivative. 

• Determined the Laplace transform of the initialization function 

• Generalized the Laplace transform of the differintegral (eq. (6.7. 1 )). 

• Applied the above theoretical results to a wide variety of applications to demonstrate the 
power of the mathematics. Showed electrical computing elements that are theoretically 
capable of simulating the generalized (fractional) operators, and demonstrated some of 
the mathematical properties through use of these analog elements. 
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• Applied the theoretical results to solve applications, in instrumentation, 
electrochemistry, material creep, filtering, and viscoelastically damped vibration. 

• Proposed a variable structure (order) fractional differintegral, provided motivation for 
its development, and proposed three particular Riemann-Liouville based 
implementations. Examined these for adherence to the criteria, and determined related 
Laplace transforms. Proposed Griinwald based variable structure definitions. (The 
variable structure concept appears to strengthen the connection between the fractional 
calculus and convolution theory.) 

The issue of when to use d' or d m versus D' or D m , respectively, will be dictated by the 
application at hand and the nature of the required initialization (value of i j / ). In an analytical 

setting it is the analyst’s prerogative in relation to the task at hand. 

While the above development was done with t, inferring time as the independent variable, 
the results may be readily extended to a spatial independent variable. 
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Appendix A 

Application Of Dirichlet’s Formula 


This appendix augments the proof of criterion 5. It is basically a restatement of the work of 
Ross (1974a and 1974b). It is desired to validate 

.dr'f(') = ,<',<'/(') = dr, 

= r(^>) | \ ( ' ~ z) "'f 1 ” r ' / ( r ') rfr ' dx ■ < A1) 

Now by a special case of the Dirichlet formula applied over the triangular region shown in 
figure A- 1 , this equation can be written as 

=7^1 j('- T )">- T .r7( t ' ') dTdT " <A - 2) 

or 

= \ f{T 1 ' T) ‘" ~ T,)V " dTdT ’- <A,3) 

Now for the inner integral let 

T — T | / \ 

y= => T=(/-T,)y + r, 

t-T , 

and dr = (t - r ,)<iy 
at the limits of integration 
t = r , => v = 0 and r = t => y = 1 . 

Therefore the inner integral of equation (A. 3) becomes 

A 1-b 

= ,)' 1 dr = (/-T,) J (l — v)" ' y v 'dy. (A.4) Figure A- 1— Reference triangle. 
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Appendix B 

Laplace Transform Of Multiple Integer Order Integrals 


This appendix will derive an expression for the Laplace transform of multiple integer order 
integrals for comparison with fractional transforms. Then, the problem is to determine 


L{g(r)} = L J j J J /(/„) dt n dt n _ t dt 2 dt t . 


(B.l) 


For notational convenience let 


J 

g ( ( j) = )<*;+! ’ for j = (B.2) 


and 


$('„-.) = }f(h) d h- 

a 

Then starting from the outside, 

t{«(')}= ijj *(>,)*,! 

■ 7 )}+;!*(',)*. -7^',)}+ 

a 

Now replacing g(/,) using equation (B.2) yields 


Ms(0}=- 


L |j s(h) dl 


S 


C l ' J S S 


Repeating the process a total of n times yields the final result, 

t h 'mi i | 


(B.3) 


(B.4) 


(B.5) 


(B.6) 

(B.7) 

(B.8) 


j j j f{t„) dt n dt «- 1 d h dt \ [= ^T L {/(0} + Zt7 » = 1,2,3,--- * (B.9) 

l a a a a 1 f = I 


where the coefficients are given by 


u 

c , =/#(';)<*. 


(B. 10) 
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Appendix C 

Variable Structure (Order) Differintegral Griinwald Basis 


A variable structure differintegral may also be formulated based on the Griinwald definition, 
equation (4.2. 1). Let AT = (/ - a) / N in equation (4.2. 1 ), and limit consideration to q => q(t) or 
q => q[t — jAT ) . Then, expressing q generally as q => q(t, jAT ) , a generalized Griinwald form 
may be written as 


Df ,) f{t)= lim XAT^ T) 


r(- q r ,(i,j&r))r(j+i) 




(Cl) 


It is observed that q(t, 7'A7’)occurs three times in this expression, in the exponent q E , in the 
numerator q N , and in the denominator q D . Thus the number of permutations possible with this 
form is eight, which is comparable with the nine possible permutations using the Riemann- 
Liouville form for the differintegral. Table C-l shows the possible permutations with the 
Griinwald formulation. 


Table C-l Permutations of 9 


Permutation 

l 

2 

3 

4 

5 

6 

1 

8 

9 , (/. /AT ) = 

t 

t 

t 

t 

t - jAT 

t - jAT 

t - jAT 

t - jAT 

<?*(/- i A f)= 

t 

t 

t - jAT 

t- jAT 

t - jAT 

t- jAT 

t 

t 

II 

<1 

t 

t - jAT 

t 

t - jAT 

t - jAT 

t 

t - jAT 

t 


Consider two Griinwald based (eq. (C.l)) evaluations of u Df ,] f(t ) , for a common q(t ) , 
with one evaluation for a < t < t { and another for a < t <t 2 where t 2 > /, . In particular, consider 
q(t)as shown in figure C-l (a), where q(t) steps from constant value </(/,) to the constant value 
q(t 2 ) at time t = . The evaluation to t = 1 2 can be viewed as an evaluation from t = a to , 

summed to an evaluation from / = /, to /, . In view of the hereditary nature of the fractional 
differintegral, the evaluation to / = /, is part of the history in the evaluation a to t 2 ; thus it is 


apparent that for the evaluation of a Df' ] f{t) 


to t = /-, 


based on the value of q(t, jAT ) = q(t) (i.e., 


all permutations except permutation 5, table C-l) will yield an undesirable result. This is so 
because, over the / < /, part of the evaluation the value of q(t, jAT) - q(t) - q(t 2 )wi\\ be used in 
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the summation, essentially changing the history (figure C-l(b)). Thus it appears that the only 
(potentially) satisfactory permutation from this set of possibilities (table C-l) is permutation 5, 
namely. 


N - 1 


a Df' ) f{t)= limXAT 


vAt-fiT) 


N-* oo 

AT-^O J~ 


rU-qA'-jAT)) 
r(-q„(t-jAT))r(j + \) 


f('-jAT). 


(C.2) 


This definition parallels equation (8. 1.4.4), and adds to the credibility of the convolution related 
variable structure form. 

i 

The intuitive nature of the Grunwald form, in terms of the available conceptualizations 
allows the hope that a form of it may be evolved which would satisfy composition. Detailed 
studies of the Grunwald permutations have not yet been performed. 




• •• 43210 = j 


(b) Functions versus j 


Figure C-l . — Variable q differintegration. (a) Functions versus t. (b) Functions versus j. 


NASA/TP— 1998-208415 


102 




References 


Apostol, T.M.: Mathematical Analysis: A Modem Approach to Advanced Calculus. 
Addison-Wesley Publishing Co., 1957. Library of Congress Card Number 57-8707. 

Bagley, R.L.: The Thermorheologically Complex Material. Int. J. Engng. Sci., vol. 29, no. 7, 
1991, pp. 797-806. 

Bard, A.J.; and Faulkner, L.R.: Electrochemical Methods. John Wiley and Sons, New York, NY, 
1980. ISBN 0-47-05542-5. 

Bland, D.R.: The Theory of Linear Viscoelasticity. Pergamon Press, Oxford, England, 1960. 
LCCCN-59- 14489. 

Churchill, RV.: Operational Mathematics. Second ed., McGraw Hill, New York, NY, 1958. 
Library of Congress Card Number 58-6678. 

Davis, H.T.: The Theory of Linear Operators: From the Standpoint of Differential Equations of 
Infinite Order. The Principia Press, Bloomington, IN, 1936. 

Gabel, R.A.; and Roberts, R.A.: Signals and Linear Systems. John Wiley and Sons, New York, 
NY, 1973. ISBN 0-471-28900-0. 

Hartley, T.T.; and Lorenzo, C.F.: Insights into the Initialization of Fractional Order Operators 
Via Semi-Infinite Lines. NASA TM-208407, 1998. 

Hartley, T.T.; Lorenzo, C.F.; and Qammer, H.K.: Chaos in a Fractional Order Chua’s System. 
IEEE Trans. Circuits Sys. I: Fundamental Theory and Applications, vol. 42, no. 8, Aug. 1995, 
pp. 485—490. 

Kalia, R.N., ed.: Recent Advances in Fractional Calculus (Global Research Notes Series in 
Mathematics). Global Publishing Co., Sauk Rapids, MN, 1993. ISBN 0-9638155-0-4. 

Kiryakova, V.S.: Generalized Fractional Calculus and Applications (Pittman Research Notes in 
Mathematics Series, 301). Longman Scientific and Technical, Essex, England and John Wiley 
and Sons, New York, NY, 1994. ISBN 0-582-219779, ISSN 0269-3674. 

McBride, A.C.: Fractional Calculus and Integral Transforms of Generalized Functions. 

Pitman Advanced Publishing, San Francisco, CA, 1979. 

McBride, A.C.; and Roach, G.F.: Fractional Calculus. Pittman Publishing Limited 
(Research Notes in Mathematics Series, 138). Longman Scientific and Technical, 

Essex, England and John Wiley and Sons, New York, NY, 1985. ISBN 0-273-08753-3. 

Mikusinski, J.: Operational Calculus. Pergamon Press, New York, NY, 1959. 

Miller, K.S.; and Ross, B.: An Introduction to the Fractional Calculus and Fractional Differential 
Equations. John Wiley and Sons, New York, NY, 1993. ISBN 0-471-58884-9. 


N AS A/TP— 1998-208415 


103 



Nishimoto, K.: An Essence of Nishimoto’s Fractional Calculus — Integrations and 
Differentiations of Arbitrary Order. Descartes Press Co., Koriyama, Japan, 1991. 

Oldham, K.B.; and Spanier, J.: The Fractional Calculus — Theory and Applications of 
Differentiation and Integration to Arbitrary Order. Vol. 1 1 1 Mathematics in Science and 
Engineering, Academic Press, Inc., New York, NY, 1974. ISBN 0-12-525550-0. 

Ross, B.: The Development, Theory, and Applications of the Gamma-Function and a Profile of 
Fractional Calculus. Ph.D. Thesis, New York Univ., 1974, S94—10032. 

Ross, B.: A Brief History and Exposition of the Fundamental Theory of Fractional Calculus. 
Lecture Notes in Mathematics, vol. 457, ed., Dold, A. and Eckmann, B.: Fractional Calculus 
and Its Applications. Proceedings of the International Conference Held at the University of 
New Haven, June 1974b, ed. Ross, B. Springer-Verlag, Berlin, 1974, ISBN 3-540-071 6 1-X, 
ISBN 387-071 6 1-X. 

Samko, S.G.; Kilbas, A. A.; and Marichev, O.I.: Fractional Integrals and Derivatives — 

Theory and Applications. Gordon and Breach Science Publishers, Philadelphia, PA, 1993. 

ISBN 2-88124-864-0. 

Smit, W.; and deVries, H.: Rheological Models Containing Fractional Derivatives. Rheol. Acta, 
vol. 9, no. 4, 1970, pp. 525-534. 

Spanier, J.; and Oldham, K.B.: An Atlas of Functions. Hemisphere Publishing Corp., 
Washington, DC, 1987. ISBN 0-891 16-573-8. 

Wylie, C.R.: Advanced Engineering Mathematics. Fourth ed., McGraw Hill, New York, NY, 
1975. ISBN 0-07-072180-7. 


NAS A/TP— 1 998-2084 1 5 


104 



Application References 


Acoustics, Boundary Layer, and Turbulence 


Keller, J.J.: Propagation of Simple Non-Linear Waves in Gas Filled Tubes With Friction, 
Journal of Applied Mathematics and Physics (ZAMP), vol. 32, Mar. 1981, pp. 170-181. 

Ruelle, D.; and Takens F.: On the Nature of Turbulence. Comm. Math. Phys., vol. 20, no. 3, 
1971, pp. 167-192. 

Sugimoto, N.: ‘Generalized’ Burgers Equations and Fractional Calculus. Nonlinear Wave 
Motion. Longman Scientific and Technical, Harlow, England, 1989, pp. 162-179. 

Sugimoto, N.: Evolution of Nonlinear Acoustic Waves in a Gas-filled Pipe. Frontiers of 
Nonlinear Acoustics: Proceedings of 12 Ih 1SNA, ed., Hamilton, M.F. and Blackstock. 

Elsevier Science Publishers Ltd., London, England, 1990. 

Sugimoto, N.: Burgers Equation With a Fractional Derivative: Hereditary Effects on Nonlinear 
Acoustic Waves. J. Fluid Mech., vol. 225, 1991, pp. 631—653. 

Controls and Dynamics 

Axtell, M.; and Bise, M.E.: Fractional Calculus Applications in Control Systems. Proceedings 
of the IEEE National Aerospace and Electronics Conference, 1990, pp. 563-566. 

Bagley, R.L.: “The Initial Value Problem for Ordinary Differential Equations with Constant 
Coefficients”, Air Force Institute of Technology, AD-203 843, N89-20764, 

AFIT-TR-EN 88-1, September 1988. 

Bagley, R.L. and Calico, R.A.: “Fractional Order State Equations for the Control of 
Viscoelastically Damped Structures", Journal of Guidance, Control, and Dynamics, Vol. 14, 
No. 2, pp. 304-3 1 1 , Mar.-Apr. 1991. 

Mbodje, B.; and Montseny, G.: Boundary Fractional Derivative Control of the Wave Equation. 
IEEE Trans. Automat. Cont., vol. 40, no. 2, Feb. 1995. 

Skaar, S.B.; Michel, A.N.; and Miller, R.K.: Stability of Viscoelastic Control Systems. 

IEEE Trans. Automat. Cont., vol. 33, no. 4, Apr. 1988. 

Diffusion 

Metzler, R.; Glockle, W.G.; and Nonnenmacher, T.F.: Fractional Model Equation for 
Anomalous Diffusion. Physica A, vol. 21 1, no. 1, 1994, pp. 13-24. 

Scher, H.; and Montroll, E.W.: Anomalous Transit-Time Dispersion in Amorphous Solids. 
Physics Rev. B, vol. 12, no. 6, Sept. 1975, pp. 2455-2477. 


NASA/TP— 1998-208415 


105 



Schnieder, W.R.; and Wyss, W.: Fractional Diffusion and Wave Equations. J. Math. Phys., 
vol. 30, no. l,Jan. 1989, pp. 134-144. 

Wyss, W.: The Fractional Diffusion Equation. J. Math. Phys., vol. 27, no. 11, Nov. 1986. 

\/f p Electrical Circuit Noise 

Bittel, H.: Noise in Magnetic and Dielectric Materials. Physica B & C, vol. 83B & C, 

1976, pp. 6-13. 

Hooge, F.N.: Discussion of Recent Experiments on 1 If Noise. Physica, vol. 60, no. 1, 

1972, pp. 130-144. 

Hooge, F.N.: 1 If Noise. Physica B & C, vol. 83B & C, May 1976, pp. 14-23. 

Keshner, M.S.: 1 If Noise. Proc. IEEE, vol. 70, no. 3, Mar. 1982, pp. 212-218. 


Electrical and Electrochemistry 

Bard, A.J.; and Faulkner, L.R.: Electrochemical Methods: Fundamentals and Applications. 
John Wiley and Sons, New York, NY, 1980. ISBN 0-471-05542-5. 

Halijak, C.A.: An RC Impedance Approximant to {VsY 1 . Correspondence, IEEE Trans. Circuit 
Theory, Dec. 1964, pp. 494—495. 

Ichise, M.; Nagayanagi, Y.; and Kojima, T.: An Analog Simulation of Non-Integer Order 
Transfer Functions for Analysis of Electrode Processes. J. Electroanal. Chem. Interfacial 
Electrochem., vol. 33, no. 2, Dec. 1971, pp. 253-265. 

Kanoh, H.; Feng, Q.; Hirotsu, T.; and Ooi, K.: A.C. Impedance Analysis for Li+ Insertion of a 
Pt / A -MnO: Electrode in an Aqueous Phase. J. Electrochem. Soc., vol. 143, no. 8, Aug. 1996, 
pp. 2610-2615. 

Oldham, K.B.: Semi-integral Electroanalysis: Analog Implementation. Analytical Chemistry, 
vol. 45, no. 1, Jan. 1973, pp. 39—47. 

Sandi, G.; Winans, R.E.; and Carrado, K.A.: New Carbon Electrodes for Secondary Lithium 
Batteries. J. Electrochem. Soc., vol. 143, no. 5, May, 1996, pp. 195-198. 

Yu, Q.; and Wing, O.: Computational Models of Transmission Lines With Skin Effects and 
Dielectric Loss. IEEE Trans. Circuits Sys.-l: Fundamental Theory and Applications, vol. 41, 
no. 2, Feb. 1994, pp. 107-119. 


NASA/TP — 1998-2084 1 5 


106 



